(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organizatioii 
International Bureau 

(43) International Publication Date 
13 March 2003 (13.03,2003) 




PCT 



lillilllllllllli 

(10) International Publication Number 

WO 03/020850 A2 



(51) International Patent Classification^: ClOG 25/03, 
25/05 

(21) International Application Number: PCT/US02/28040 

(22) International Filing Date: 

4 September 2002 (04.09.2002) 



(25) Filing Language: 

(26) Publication Language: 



English 
English 



(30) Priority Data: 

60/317,158 4 September 2001 (04.09.2001) US 

(71) Applicant (for all designated States except US): THE 
REGENTS OF THE UNIVERSITY OF MICHIGAN 
[US/US]; 3003 South State Street, Wolverine Tower, Room 
2071, Ann Aibor, Ml 48109-1280 (US). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): YANG, Ralph, X 
[US/US]; 2627 Pin Oak Drive, Ann Arbor, MI 48103 (US). 



YANG, Frances, H. [US/US]; 2627 Pin Oak Drive, Ann 
Arbor, Mi 48103 (US). TAKAHASHI, Akira [JPAJS]; 
1653 Mclntyre Drive, Ann Arbor, MI 48105 (US). HER- 
NANDEZ-MALDONADO, Arturo, J. [US/US], 1755 
Mclntyre Street, Northwood 4, Ann Ari>or, MI 48105 
(US). 

(74) Agent: DIERKER, Julia, Church; Dierker & Glass- 
meyer, RC, 3331 W. Big Beaver Road, Suite 109, Troy, 
MI 48084-2813 (US). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, H, GB, GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT. LU, LV, MA, MD, MG, MK, MN. MW, 
MX, MZ, NO, NZ, OM, PH, PL, PT, RO. RU, SD, SE, SG, 
SI, SK, SL, TJ. TM, TN, TR. TT. TZ, UA, UG, US, UZ, 
VN, YU, ZA, ZM, ZW. 

(84) Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ, SD, SL. SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ. TM). 

f Continued on next page]. 



(54) Title: SELECTIVE SORBENTS FOR PURIRCATION OF HYDROCARBONS 



IT) 
00 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 




_ 



BENZENE ON ^ 
Na-Y N^/ 



THIOPHENEON 
Na-Y 



>i I iMtil I I iimil I 1 1 imil 1 t mini « « 



IlL 



1.E-06 1.E-05 1.E-04 1.E-03 1.E-02 
PARTIAL PRESSURE (atm) 



1.E-01 1.E+00 



(57) Abstract: Cu Y and Ag Y zeolites as selective sorbents for desulfurization of liquid fuels. Thiophene and benzene were used 
as the model system, and vapor phase isotherms were measured. Compared with NaY, CuY and AgY adsorbed significantly laiger 
^2 amounts of both thiphene and benzene at low pressures. It is hypothesized that this is due to n-complexation with Cu* and Ag*. On 
a per-cation basis, more thiophene was adsorbed by Cu* than by Ag*, eg. 0.92 molecule/Cu* versus 0.42 molecule/Ag* at 2x10-^ 
atm and 120 °C. Molecular orbital calculations confirmed the relative strengths of TC-complexation: thiophene>benzene and Cu*> 
Ag*. The experimental heats of adsorption for 7C-complexation are in qualitative agreement with theorical predictions. The invention 
further comprises a process and sorbents for removal of aromatics from hydrocarbons. 



wo 03/020850 A2 lilllllilllllllllllllllllllllllllilllllllllllllli 



European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, 
ES, FI, FR, GB, GR, IE, IT, LU, MC, NL, PT, SE, SK, 
TR), OAPI patent (BF, BJ, CF, CG, CI, CM, GA. GN, GQ, 
GW. ML, MR, NE, SN, TD, TG). 



For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



Published: 

— without international search report and to be republished 
upon receipt of that report 



wo 03/020850 



PCT/US02/28040 



Our Reference: UMJ-116-B(UM-2172) PATENT 
SKT RrnVE SOKB ENTS FORPITRTFICATION OF HYDROCARBONS 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional Patent application 

Serial No. 60/317,158, filed September 4, 2001. 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 
This invention was made in the course of research partially supported by a 

grant from the National Science Foundation (NSF) (Grant No. CTS-98 19008 and Grant 

No.CTS-01 381 90); and by a grant firom the Department of Energy (DOE) (Fuel Cell 

Grant No. DE-FC04-02AL67630). The U.S. government has certain rights in the 

invention. 

BACKGROUND OF THE INVENTION 
The present invention relates generally to processes for purification of 
hydrocarbons, and more particularly to adsorption processes using sorbents selective to 
sulfiir compounds and to aromatic compounds. 

Petroleum is an extremely complex mixture and consists predominantly of 
hydrocarbons, as well as compounds containing nitrogen, oxygen, and sulfur. Most 
petroleums also contain minor amounts of nickel and vanadium. The chemical and 
physical properties of petroleum vary considerably because of the variations in 
composition. 

The ultimate analysis (elemental composition) of petroleum tends to vary 
over relatively narrow limits — carbon: 83.0 to 87.0 percent; hydrogen: 10.0 to 14.0 
percent; nitrogen: 0.1 to 1.5 percent; oxygen: 0.1 to 1.5 percent; sulfur: 0.1 to 5.0 percent 
metals (nickel plus vanadium): 10 to 500 ppm. 

Crude oils are seldom used as fiiel because they are more valuable when 
refined to petroleum products. Distillation separates the crude oil into fi-actions 
equivalent in boiling range to gasoline, kerosine, gas oil, lubricating oil, and residual 
Thermal or catalytic cracking is used to convert kerosine, gas oil, or residual to gasoline, 
lower-boiling fiactions, and a residual coke. Petrochemical intermediates such as 
ethylene and propylene are primarily produced by the thermal cmcking of light 
hydrocarbon feedstocks in the presence of steam. Catalytic reforming, isomerization, 
allfylation, polymerization, hydrogenation, and combinations of these catalytic processes 
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are used to upgrade the various refinery intermediates into improved gasoline stocks or 
distillates. The major finished products are usually blends of a number of stocks, plus 
additives. 

Gasoline is a complex mixture of hydrocarbons that distills within the 
range 100 to 400^F. Commercial gasolines are blends of straight-run, cracked, reformed, 
and natural gasolines. Straight-run gasoline is recovered jfrom crude petroleum by 
distillation and contains a large proportion of normal hydrocarbons of the parafiBn series. 
Cracked gasoline is manufectured by heating cmde-petroleum distillation fiuctions or 
residues under pressure, or by heating with or without pressure in the presence of a 
catalyst Heavier hydrocarbons are broken into smaller molecules, some of vdiich distill 
in the gasoline range. Reformed gasoline is made by passing gasoline fiactions over 
catalysts in such a manner that low-octane-number hydrocarbons are molecularly 
rearranged to high-octane-number components. Many of the catalysts use platinum and 
other metals deposited on a silica and/or alumina support. Natural gasoline is obtained 
fi*om natural gas by liquefying those constituents which boil in the gasoline range either 
by compression and cooling or by absorption in oil. 

Removal of the sulfur-containing compounds is an important operation in 
petroleum refining, and is achieved by catalytic processes at elevated temperatures and 
pressures. See, Farrauto, R. J.; Bartholomew, C. H. Fundamentals of Industrial Catalytic 
Processes, Chapman and Hall, New York, 1997. The hydrodesulfiirization (HDS) 
process is eflBcient in removing thiols and sulfides, but much less eflfective for 
heterocyclic diunsaturated sulfiir compounds, such as thiophenes and thiophene 
compounds/derivatives (eg. benzothiophene and dibenzothiophene). 

New legislation will require substantial reductions in the sulfiir content of 
transportation fiiels. For example, the new U.S. Environmental Protection Agency (EPA) 
sulfiu- standards require that the sulfiir content in gasoline and diesel fiiels for on-board 
transportation will be 30 ppm and 1 5 ppm, respectively, drastically decreased from the 
current acceptable levels of several hundred ppm. Faced with the severely high costs of 
compliance, a surprising number of refiners are seriously considering reducing or 
eliminating production of on-board fiiels. See Parkinson, G., 'Diesel Desulfiirization 
Puts Refiners in a Quandary," Chemical Engineerings 2001, February issue, 37. 

During the last decade, there have been several published accounts on 
using adsorption for liquid fiiel desulfiirization. Commercially available sorbents (i.e., 
zeolites, activated carbon and activated alumina) were used in all of these studies. 
Weitkamp et al. reported that thiophene adsorbed more selectively than benzene on ZSM- 
5 zeolite. See Weitkamp, J.; Schwark, M.; Ernest, S. "^Removal of Thiophene Impurities 
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from Benzene by Selective Adsorption in Zeolite ZSM-5," J. Chem, Soc. Chem. 
Commun,, 1991, 1 1 33 . Without being bound to any theory, it is believed that this is 
because thiophene (C4H4S, also known as thiofuran) has a higher dipole moment (0.55 
debye) than benzene (non-polar), although their polarizabilities are similar. Based on this 
study. King et al studied selective adsorption of thiophene, methyl- and dimelfayl- 
thiophenes (all with one ring) over toluene and p-xylene, also using ZSM-5. See King, 
D. L; Faz, C; Flynn, T. "Desulfurization of Gasoline Feedstocks for Application in Fuel 
Reforming," SAE Paper 2000-01-0002, Soc. Automotive Engmeers, Detroit, 2000. They 
showed that thiophene was more selectively adsorbed, both based on fixed bed 
breakthrough experiments. However, the capacities for thiophene were xmfortunately 
quite low (only 1 -2% wt adsorbed at 1% thiophene concentration). Both vapor phase 
and liquid phase breakthrough experiments were done in these studies, and the results 
from two phases were consistent 

The pore dimensions of ZSM-5 are 5.2-5.6 A, Hence, organic sulfur 
compounds with more than one ring will be sterically hindered or excluded. Zeolites 
with larger pores, as well as larger pore volumes, would appear to be more desirable than 
ZSM-5 as the selective sorbents. Lideed, results of Salem and Hamid indicated that 13X 
zeolite ais well as activated carbon had much higher sorption capacities for sulfur 
compounds. See Salem, A. S. H.; Hamid, H. S. "Removal of Sulfur Compounds from 
Naphtha Solutions Using Solid Adsorbents," Chem. Eng. Tech., 1997, 20, 342. Based on 
the data of Salem and H^mid, Id., the capacity for sulfur compounds by 13X zeolite was 
approximately an order of magnitude higher than that of ZSM-5, v/hen compared with the 
data of King et al. (cited above) extrapolated to the same conditions. 

Activated alumina (Alcoa Selexsorb) has been used in an adsorption 
process by Irvine. See, for example, U. S. Patent No. 5,730,860 issued to R.L. Irvine in 
1998, entitled "Process for Desulfurizing Gasoline and Hydrocarbon Feedstocks." 

No direct comparison has been made among these commercial sorbents. 
Their experiments were mostly done in fixed bed adsorbers, by measuring the 
breakthrough capacities. Based on the literature, the large pore zeolites (NaX or NaY) 
are about the same as activated carbon and alumina, in terms of adsorption of thiophene. 

Essentially all industrial adsorption processes are based on van der Waals 
interactions between the sorbate and the sorbent Chemical bonds have yet to be 
exploited. 

Thus, it would be desirable to provide an adsorption process for 
selectively removing sulfur compounds from liquid fuels at ambient temperature and 
pressure, thereby advantageously leading to a major advance in petroleum refining. It 
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would further be desirable to provide highly selective sorbents for this process, thereby 
overcoming the drawbacks of current commercial sorbents, which are not desirable for 
this application. 

SUMMARY OF THE INVENTION 

The present invention addresses and solves the above-mentioned 
drawbacks, and substantially meets the desiderata and advantages enumerated above, by 
providing a process for removing thiophene and thiophene compounds fix)m liquid &el. 
The method comprises the step of contacting the liquid fuel with an adsorbent which 
preferentially adsorbs the thiophene and thiophene compounds, at a selected temperature 
and pressure, thereby producing a non-adsorbed component and a thiophene/thiophene 
compound-rich adsorbed component. The adsorbent comprises any ion-exchanged 
zeolite, but in a preferred embodiment, the zeolite is selected from the group consisting of 
zeolite X, zeolite Y, zeolite LSX, and mixtures thereof The zeolite has exchangeable 
cationic sites, and at least some of the sites have d-block transition metal cation present 
The preferential adsorption occurs by Ti-complexation. 

In a further embodiment, the invention comprises a process for removing 
thiophene and thiophene compounds from liquid fiiel. The method comprises the step of 
contacting the liquid fiiel with an adsorbent which preferentially adsorbs the thiophene 
and thiophene compounds, at a selected temperature and pressure, thereby producing a 
non-adsorbed component and a thiophene/thiophene compound-rich adsorbed 
component. The adsorbent comprises a canier having a sur&ce area, the carrier having a 
d-block transition metal compound dispersed on at least some of the sur&ce area. The 
metal compound releasably retains the thiophenes. The carrier comprises a plurality of 
pores having a pore size greater than the effective molecular diameter of the 
thiophene/thiophene compounds. The method may ftirther comprise the step of changing 
at least one of the pressure and temperature to thereby release the thiophene/thiophene 
compound-rich component from the adsorbent 

In a fiirther embodiment, the invention comprises a process for selective 
removal of aromatics from hydrocarbons. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Objects, features and advantages of the present invention will become 
apparent by reference to the following detailed description and drawings, in v^ich: 

Fig. 1 is a graph depicting equilibrium isotherms of benzene and thiophene 
on NaY (Si/Al = 2.43) at 120^C; with the lines being fittings by Langmuir-Freundlich 
isotherm equation; 
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Fig. 2 is a graph depicting equilibrium isothenns on AgY (56 AgVunit 
cell) and CuY (14 Cu^nit cell) at 120^C; 

Fig. 3 is a graph depicting thiophene/benzene and thiophene/n-octane 
breakthrough in AgY, NaY zeolite at room temperature; 

Fig. 4 is a graph depicting cumulative pore volumes of sorbents used in 

Experiment C; 

Fig. S is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on Na-Y (Si/Al=2.43) at 120 "^C and 1 80 ""C; curves are fitted 
with Dubinin-Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 

Fig. 6 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on Ag-Y (Si/Al=2.43) at 1 20**C and 1 80°C; curves are fitted with 
Dubinin-Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; fitted 
curves are not shown for benzene adsorption at 1 80°C because the artificial crossovers to 
the curves for thiophene at 1 80°C are observed; 

Fig. 7 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on Cu-Y (Si/Al=2.43) at 90''C and 120^*0; curves are fitted with 
Dubinin-Astakhov (solid line) and Langmuir-Freundlich (dotted line) isothenns; 

Fig. 8 is a graph depicting pure component equilibrium isotherms of 
thiophene on Ag-X (Si/Al=1.25) and Ag-Y(Si/Al=2.43) at 120**C; 

Fig. 9 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on H-USY (Si/Ai=195) at 90°C and 120*^0; curves are fitted with 
Dubinin-Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 

Fig. 10 is a graph depicting pure component equilibrium isotherms of 
benzene and tiiiophene on Na-ZSM-5 (Si/Al=10) at 90°C and 120*='C; curves are fitted 
with Dubinin-Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 

Fig. 1 1 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on activated carbon (Type PCB) at 90°C and 120°C; curves are 
fitted with Dubinin-Astakhov (solid line) and Langmuir-Freundlich (dotted line) 
isotherms; 

Fig. 12 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on modified activated alumina (Selexsorb CDX) at 90°C and 
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120^C; curves are fitted with Dubinin-Astakhov (solid line) and Langmuir-Freundlich 
(dotted line) isothenns; 

Fig. 13 is a graph comparing equilibrium adsorption isothenns of thiophene 

at 120^C; 

Figs. 14-21 are graphs depicting various breakthrough curves; 

Fig. 22 is a graph depicting pure component equilibrium isothenns for 
benzene and cyclohexane adsorption on H-USY at 120 and 1 80 ^C; curves are fitted with 
Dubinin-Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isotherms; 

Fig. 23 is a graph depicting pure component equilibrium isotherms for 
benzene and cyclohexane adsorption on Na-Y at 120 and 180 °C; curves are fitted with 
Dubinin-Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isothenns; 

Fig. 24 is a graph depicting pure component equilibriiun isotherms for 
benzene and cycloheTcane adsorption on Ag-Y at 1 20 and 1 80 ^C; curves are fitted with 
Dubinin-Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isotherms; 

Fig. 25 is a graph depicting pure component equilibrium isotherms for 
benzene and cyclohexane adsorption on Pd-Y at 120 °C; curves are fitted with Dubinin- 
Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isotherms; 

Fig. 26 is a graph depicting pore size distributions of Na-Y, Ag-Y, Pd-Y 

and H-USY; 

Fig. 27 is a graph comparing experimental and GCMC simulation for 
adsorption of benzene on Na-Y; 

Fig. 28 is a graph comparing experimental and GCMC simulation for 
adsorption of benzene on H-USY; and 

Fig. 29 is a graph comparing experimental and GCMC simulation for 
adsorption of benzene on Ag-Y. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention is predicated upon the unexpected and fortuitous 
discovery that thiophene and thiophene compounds are adsorbed slightly more selectively 
via n -complexation than is benzene. This is quite counter-intuitive, as it would be 
expected that benzene, having more double bonds (3) and more k electrons than 
thiophene (2 double bonds), v^ould be more selectively adsorbed viaTc-complexation. An 
example of a compound having more double bonds being more selectively adsorbed than 
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a compound having fewer double bonds may be found in U.S. Patent No. 6,21 5,037, 
issued to Padin, Munson and Yang in 2001 and entitled, "Method for Selective 
Adsorption of Dienes." 

Without being bound to any theory, it is believed that this counter- 
intuitive, slightly higher selectivity for thiophene/thiophene derivatives/compounds may 
be explained by the following theory. The sulfur atom in the five-member thiophene ring 
has more electrons than the carbons. As such, the S atom, with its available electrons and 
relatively strong attraction, may be aiding in tiie tc -complexation bonding, thus 
contributing to the higher selectivity of the present sorbents for thiophenes over benzene. 
How much the sulfur could contribute to n -complexation bonding, however, is not 
predictable. 

The present invention is a startling and usefiil discovery, in that the 
sorbents and process of the present invention adsorb more than 1 mmol/gram of 
thiophene at 10'^ atm. The present inventive process utilizing sorbents which adsorb via 
7C-complexation may be used to remove sulfur compounds (for example, heterocyclic 
diunsaturated sulfur compounds) from hydrocarbons such as liquid fuels (including, but 
not limited to gasoline, diesel fuels, coal and shale derived liquid fuels, methanol, and the 
like). Liquid fiiels have a high concentration of aromatics, such as eg. benzene, and a low 
concentration of thiophene and thiophene compounds (in the parts-per-million (ppm)). 
Thus, it is one important aspect of the present invention that the inventive process 
adsorbs the thiophenes at least slightly more selectively than benzene. 

This is a substantial advance over currently available commercial sorbents, 
which adsorb little, if any, thiophene or thiophene compounds at low concentrations. For 
example, activated carbon was recently tested by the present inventors, and was found to 
behave similarly to Na-Y (discussed in fiirther detail below), in that it adsorbed very little 
(below 1/10 mumol/gram) benzene and thiophene at 10"^ atm. 

Without being bound to any theory, it is believed that the higher sorbent 
capacity of embodiments of the present invention may be due to a sorbent pretreatment 
method of an embodiment of the present invention wherein the sorbent is calcined at 
350°C and cooled in dry atmosphere. 

The present inventive process may advantageously be run at ambient 
temperature and pressure, which is highly desirable for a variety of reasons. It is much 
less energy consuminjg to run processes at ambient temperature and pressure. Further, 
conventional catalytic desulfurizing processes, processes are typically run at high 
temperatures, eg. 700**C, and high pressures, causing the (very expensive) catalyst to be 
continually deactivated. Further, the reactor size of a conventional process would need to 
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be exponentially increased over present reactor sizes in order to remove thiophenes lower 
than the currently acceptable several hundred ppm. Still further, the present inventive 
process may be used as a first line desulfurizing process, and/or as a clean up 
desulfurizing process to remove sulfur compounds missed by conventional processes, 
such as those processes discussed above (eg. hydrodesulfiirization). It could be veiy 
useful to have the option to add the present invention as a downstream "clean up*' 
desulfurizing process, in that it may not be necessary to revamp current refining 
processes upstream fi*om the present inventive process. This could solve the cost- 
prohibitive problems encoimtered by some refiners. See Parkinson, G«, (cited above). 

The sorbents of the present invention, one example of which is Cu(I)Y, 
may be substantially fiilly regenerated by first air oxidation (e.g., at about 350**C), 
followed by auto-reduction in an inert atmosphere (e.g., at about 450*^0). 

The process and sorbents of the present invention may have quite 
significant cost benefits on the future energy picture. For example, for a fiill tank of gas, 
it is estimated that it would cost only about 70 cents to remove the sulfur firom the current 
350 ppm to below 5 ppm. If some blending were done, to meet the new (2006) EPA 
standard, then it is estimated that it would cost only about 12 cents. These costs are 
based on the assumption that the sorbent is disposable - no regeneration. However, it is 
estimated that the cost could be considerably less if the sorbent were regenerated (for 
example, as described above, ie. air oxidation followed by auto-reduction in inert 
atmosphere). 

Although the process of the present invention has specifically tested Cu-Y 
and Ag-Y (as well as other sorbents), it is to be understood that Type X zeolites may in 
some cases be as good as, or better zeolites than Y zeolites, since more cations are 
available in X zeolites. Further, it is to be understood that other zeolites are contemplated 
as being within the scope of the present invention. Still j&irther, it is to be understood that 
all d-block transition metals and/or their ions may be used in place of the copper or silver, 
as it is believed that any d-block element will form n -complexation bonds with thiophene 
and thiophene compounds. In particular, it is believed that Ni^^ and Pd^^ would be as 
effective as Cu"*" and Ag^ Further, the metals do not need to be ion-exchanged, but rather 
may be dispersed (monolayer dispersion, island dispersion, etc.) on a carrier (such as, for 
example, silica, alumina, etc.) by any suitable method. 

In an embodiment of the present invention, the process for removing 
thiophene and thiophene compounds fi-om liquid fuel comprises the step of contacting the 
liquid fuel with an adsorbent which preferentially adsorbs the thiophene and thiophene 
compounds, at a selected temperature and pressure, thereby producing a non-adsorbed 
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component and a thiophene/thiophene compound-rich adsorbed component It is to be 
understood that in the "thiophene/thiophene compound-rich adsorbed component" term 
as used herein, the thiophenes may not be the greatest amount of adsorbate. For example, 
in the case of liquid fuels having a very high concentration of aromatics, and a very low 
concentration of thiophenes, the adsorbed phase will actually contain more benzene than 
thiophene, because there is much more benzene than thiophene in the solution. The 
adsorbed phase may have thiophene/benzene at a ratio of about 1/2. However, that would 
be enough to advantageously purify the solution* 

The adsorbent comprises any ion-exchanged zeolite, but in a preferred 
embodiment, the zeolite is selected from Ihe group consisting of zeolite X, zeolite Y, 
zeolite LSX, and mixtures thereof. The zeolite has e^f changeable cationic sites, and at 
least some of the sites have d-block transition metal cation present The preferential 
adsorption occurs by 7t-complexation. 

In a further embodiment, the invention comprises a process for removing 
thiophene and thiophene compounds from liquid fuel The method comprises the step of 
contacting the liquid fiiel with an adsorbent which preferentially adsorbs the thiophene 
and thiophene compounds, at a selected temperature and pressure, thereby producing a 
non-adsorbed component and a thiophene/thiophene compound-rich adsorbed 
component The adsorbent comprises a carrier having a surface area, the carrier having a 
d-block transition metal compound dispersed on at least some of the sur&ce area. The 
metal compound releasably retains the thiophenes. The carrier comprises a plurality of 
pores having a pore size greater than the effective molecular diameter of the 
thiophene/thiophene compounds. The method may frirther comprise the step of changing 
at least one of the pressure and temperature to thereby release the thiophene/^iophene 
compound-rich component from the adsorbent. 

As suggested by King, chemical complexation bonds are generally 
stronger than van der Waals interactions (thus giving rise to higher selectivities), yet 
weak enough to be reversible (i.e., to be broken by simple engineering means). See, for 
example. King, C. J., "Separation Processes Based on Reversible Chemical 
Complexation," Handbook of Separation Process Technology^ R. W. Rousseau ed., 
Wiley, New York, 1987, Chap. 15. Therefore, a tremendous opportunity exists for 
developing new sorbents (and new applications in separation and purification) by using 
weak chemical bonds, including various forms of complexation bonds. Based on the 
principles of Tc-complexation, the inventor(s) have already developed a number of new 
sorbents for a number of applications. These include sorbents for: (a) oldSn/paraffin 
separations (see, for example, Yang, KT.; Kikkinides, E. S. "New Sorbents for Olefin- 
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Paraffin Separations by Adsorption viajr -Complexation,"^CA£ J., 1995, 41, 509; 
Rege, S. U.; Padin, J; Yang, R T. "Olefm-Paraffin Separations by Adsorption: 
Equilibrium Separation by 7t -Complexation vs. Kinetic Separation," ^C/lB J., 1998, 44, 
799; Huang, a Y.; Padin, J.; Yang, R T. ''Ab Initio Effective Core Potential Study of 
Olefin/Paraffin Separation by Adsorption via % -Complexation: Anion and Cation Effects 
on Selective Olefin Adsorption," J. Phys. Chem. B., 1999, 103, 3206; Huang, H. Y; 
Padin, J,; Yang, R T. "Comparison of 7c-Complexations of Ethylene and Carbon 
Monoxide and Cu^ and Ag^," Ind, Eng. Chem, Res., 1998, 38, 2720; and Padin, J.; Yang, 
R T. "New Sorbents for Olefin-Paraffm Separations by Adsorption via % -Complexation: 
Synthesis and Effects of Substrates," Chem, Eng. Set, 2000, 55, 2607.); (b) diene/olefin 
separation or purification (i.e., removal of trace amounts of dienes fi-om olefins) (see 
Padin, J.; Yang, R T.; Munson, C. L. "New Sorbents for Olefin-Paraffm Separations and 
Olefin Purification for C4 Hydrocarbons," Ind. Eng. Chem. Res., 1999, 38, 3614; 
Jarayaman, A; Yang, RT.; Munson, C. L.; Chinn, D. "Deactivation of tt -Complexation 
Adsorbents by Hydrogen and Rejuvenation by Oxidation," Ind. Eng. Chem. Res., 2001; 
and Takahashi, A.; Yang, RT.; Munson, CX.; Chinn, D. "Cu(I)-Y Zeolite as a Superior 
Adsorbent for Diene/Olefin Separation," Langfnuir, 2001.); and (c) aromatics/aliphatics 
separation and purification (i.e, removal of trace amounts of aiomatics fi-om aliphatics) 
(see Takahashi, A; Yang, F. H.; Yang, R T. "Aromatics/Aliphatics Separation by 
Adsorption: New Sorbents for Selective Aromatics Adsorption by n -Complexation," 
Ind. Eng. Chem. Res., 2000, 39, 3856; and Takahashi, A.; Yang, R T. "New Adsorbents 
for Purification: Selective Removal of Aromatics," yl/CAEJ:, submitted (2001)). For 
purification, vapor phase results can be applied directly to liquid phase applications. In 
the present invention, the inventors have used molecular orbital calculations to obtain a 
basic understanding for the bonding between the sorbates and sorbent surfaces, and 
further, to develop a methodology for predicting and designing n -complexation sorbents 
for targeted molecules (see, for example, Huang, H.Y.; Padin, I; Yang, R T. ''Ab Initio 
Effective Core Potential Study of Olefin/Paraffm Separation by Adsorption via 7t - 
Complexation: Anion and Cation Effects on Selective Olefin Adsorption,",/. Phys. Chem. 
B., 1999, 103, 3206; and Huang, H. Y; Padin, J.; Yang, R T. "Comparison of tc - 
Complexations of Ethylene and Carbon Monoxide and Cu"^ and Ag^," Ind. Eng. Chem. 
Res., 1998,35,2720). 

In the present invention, the inventors obtained first results on new 
sorbents for desulfurization using Ti-complexation. CuY and AgY were the sorbents, and 
thiophene and benzene were used as the model system. Since NaX (i.e., 13X) orNaY is 
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among the best commercially available sorbents for thiophene, NaY was used as the 
reference for comparison. 

The sulfur breakthrough curves with AgY and NaY zeolites (Figure 3), 
with a liquid feed of thiophene in octane, shows AgY to be an excellent sorbent. 

Further experiments were run with Cu(I)Y. Surprisingly, even though the 
results witfi AgY were very good, the results with Cu(I)Y were better than those of AgY. 
The sorbent capacity of Cu(I)Y is much higher than that of AgY, and very pure octane 
was obtained for extended periods of time. The sul&r capacity of Cu(I)Y zeolite was 
foimd to be about 2.55 mmol/g (or, 21.4 wt %). This is an extremely high capacity 
because the experiments were removing ppm levels of sulfur. 

Ag-exchanged &ujasite was claimed for thiophene removal in U.S. Patent 
No. 4,188,285 issued to Michlmayr. In sharp contrast to the present invention, the *285 
preferred temperature for adsorption was 200-3 5 0°C, and the sorbent was not dehydrated 
by heat-treatment prior to adsorption. Also in sharp contrast to the present invention, the 
*285 sulfur capacities were veiy low, at about 0.07-0.15 wt% for Ag-Y. In the *285 
patent, the highest sulfur capacity (of 0.2 wt%) was obtained with the lowest Ag content, 
with USY zeolite. The *285 sorbent was apparently intended for bonding the sulfur atom 
with Ag. It was not intended for 7C-complexation. In contrast, the sorbents of the present 
invention, for example, the Cu(I)Y sorbent has a sulfur capacity > 1 00 times larger. 

Vansant et al. in European Patent Publication No. 0 275 855 investigated 
Cu(II)Y for thiophene removal. In contrast to the present invention, the *855 Cu^^ 
exchanged Y zeolite was purposely heat-treated in air (to 200-550X) to maintain the 
Cu^^ in the divalent state. Thus, in contrast to the present invention, in the *855 
publication adsorption by 7i-complexation was not intended. The '855 publication 
discloses low thiophene capacities (the highest was 1.6 wt%). In shatp contrast, the 
thiophene capacities of the present invention were about 21 ,4 wt%. 

Figs. 14-21 are graphs depicting various breakthrough curves. 

To further illustrate the present invention, the following examples are 
given. It is to be understood that these examples are provided for illustrative purposes 
and are not to be constmed as limiting the scope of the present invention. 

EXAMPLES 

EXPERIMENT A 
Sorbent Preparation 

Cu"*'-Y was prepared by ion exchange of Na-Y zeolites (Si/Al=2.43, 56 Al 
atoms/u.c, Strem Chemical) with Cu(N03)2 followed by reduction of Cu^^ to Cu**^, since 
CuOO is not soluble in water. First, as-received Na-Y was exchanged twice using excess 
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amounts (10-fold cation-exchange-capacity (CEC) assuming that one Cu^^ compensates 
two aluminum sites) of 0.5 M Cu(N03)2 at room temperature for 24 hours. After the 
exchange, the zeolite suspension was filtered and washed with copious amount of de- 
ionized water. The product was dried at lOO^C overnight Several groups have reported 
reduction of Cu^^ to Cu"*^ in zeolite in He (i.e., auto-reduction, see Takahashi, A.; Yang, 
RT.; Munson, C.L.; Chiim, D. "Cu(I)-Y Zeolite as a Superior Adsorbent for 
Diene/Olefin Separation," Langmuir, 2001; and Larsen, S. C; Aylor, A.; Bell, A. T.; 
Reimer, J. A ''Electron Paramagnetic Resonance Studies of Copper Ion-Exchanged 
ZSMrS:' J. Phys.Chem., 1994, 98, 11533). In this study, reduction of Cu^ to Cu^ was 
carried out in He only at 450°C for 1 hour. Ag^ ion-exchanged Y-zeolite (Ag-Y) was 
prepared at room temperature for 24 h in the same manner as Cu^"*^ exchange, using 5-fold 
excess AgpSTOs (O.IM). 

Single component isotherms for benzene and thiophene were measured at 
120°C using standard gravimetric methods. A Shimadzu TGA-50 automatic recording 
microbalance was employed. Helium (Pre-purified grade, Metro welding, 99.995%) was 
used as the carrier gas and was first passed through two consecutive gas-wash bottles (to 
ensure saturation), which contained benzene (HPLC grade, Aldrich, 99.9+%) or 
tiiiophene (Aldrich, 99+%). After diluting the concentration to the desired value by 
blending witih additional helium, the mixture was directed into the microbalance. 
Molecular Orbital Computational Details 

Molecular orbital (MO) studies on the n -complexation bonding for 
benzene and sorbent surfiices had been investigated recently. See Takahashi, A.; Yang, 
F. H.; Yang, R T. "Aromatics/Aliphatics Separation by Adsorption: New Sorbents for 
Selective Aromatics Adsorption by n -Complexation," Ind. Eng. Chem. Res., 2000, 39, 
3856. In this work, similar MO studies were extended to thiophene and sorbent surfaces. 
The Gaussian 94 Program in Cerius2 molecular modeling software (see Bowie, J. E., 
Data Visualization in Molecular Science: Tools for Insight and Innovation; Addison- 
Wesley Pub. Co.: Reading, Mass., 1995; chapter 9) from Molecular Simulation, Inc. was 
used for all calculations. See Frisch, M. J.; Trucks, G.W.; Schlegel, H. B.; Gill, P. M.W.; 
Johnson, B. G.; Robs, M. A.; Cheeseman, J. R; Keith, T.; Petersson, G. A,; Montgomery, 
J. A.; Raghavachari, K.; Al-Laham, M. A.; Zakrzewski, V. G; Ortiz, J. V.; Foresman, J. 
B.; Peng, C. Y.; Ayala, P. Y.; Chen, W.; Wong, M. W.; Andres, J. L.; Replogle, E. S.; 
Gomperts, R; Martin, R L.; Fox, D. J.; Binkley, J. S.; DeFrees, D. J.; Baker, J.; Stewart, 
J. P.; Head-Gordon, M.; Gonzalez, C; Pople, J. A. Gaussian 94, Revision B.3, Gaussian, 
Inc., Pittsburgh, PA., 1995. MO calculations for thiophene and sorbent surfeces were 
performed at the Hartree-Fock (HF) and density fiinctional theory (DFT) level using 
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effective core potentials (ECPs). See Hay, P. J.; Wadt, W.R. "Ab Initio Effective Core 
Potential for Molecular Calculations. Potentials for the Transition Metals Atoms Sc to 
Hg," J. Chem. Phys. 1985, 82, 270; Wadt, W. R.; Hay, P. J. "Ab Initio Effective Core 
Potential for Molecular Calculations: Potentials for Main Group Elements Na to Bi," J. 
Chem. Phys. 1985, P2, 284; and Hay, P. J.; Wadt, W. R "Ab Initio Effective Core 
Potential for Molecular Calculations: Potentials for K to Au Including the outennost core 
orbilals," J. Chem. Phys. 1985, 52, 299, 

The LanL2DZ basis set (see Russo, T. V.; Martin, R. L.; Hay, P. J. 
"Effective Core Potentials for DFT Calculations," J. Phys. Chem. 1995, 99, 17085) is a 
double-^ basis set containing ECP representations of electrons near the nuclei for post- 
third-row atoms. The reliability of this basis set has been confirmed by the accuracy of 
calculation results as compared with experimental data. Therefore, the LanL2DZ basis 
set was employed for both geometry optimization and natural bond orbital (NBO) 
analysis. 

The restricted Hartree-Fock (RHF) theory at the LanL2DZ level basis set 
was used to determine the geometries and the bonding energies of thiophene on AgCl and 
CuCl. The simplest models with only a single metal chloride interacting with a thiophene 
molecule were chosen for n -complexation studies. The optimized structures were then 
used for bond energy calculations according to the following expression: 

Eads ~ Eadsorbate + Eadsarbent Eadsoibent-edsoifoate (1) 

where Eadsoibate is the total energy of thiophene, Eadsoibent is the total energy of the bare 
adsorbent, i.e. the metal chloride, and Eadsorbenvadsorbaie is the total energy of the 
adsorbate/adsorbent system. A higher value of Eads corresponds to a stronger adsorption. 
Natural Bond Orbital (NBO) 

The optimized structures were also used for NBO analysis at the 
B3LYP/LanL2DZ level. The B3LYP (see Becke, A D. 'Density Functional 
Thennochemistry. IE. The Role of Exact Exchanges," J. Chem. Phys. 1993, 98, 5648) 
approach is one of the most useful self-consistent hybrid (SCH) approaches (see Becke, 
A D.; "A New Mixing of Hartree-Fock and Local Density-Functional Theories," J. 
Chem. Phys. 1993, PS, 1372), it is Beck's 3-parameter nonlocal exchange fiinctional (see 
Becke, A D.; 'T)ensity-Functional Thentnochemistiy. XL TheEffect of IhePerdew-Wang 
generalized-gradient correlation correction," J. Chem. Phys. 1992, 97, 9173) with 
nonlocal correlation functional of Lee, Yang and Parr. See Lee, C; Yang, W.; Parr, R.G. 
"Development of the CoUe-Salvetti Correlation-Energy Formula into a Functional of the 
Electron Density," Phys. Rev., 1988, B37, 785. 
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The NBO analysis performs population analysis that pertains to localized 
wave-function properties. It gives a better description of the electron distribution in 
compounds of high ionic character, such as those containing metal atoms. See Reed, A. 
E.; Weinstock, RB.; Weinhold, F. "Natural Population Analysis," J. Chem. Phys, 1985, 
83 (2), 735. It is known to be sensitive for calculating localized weak interactions, such 
as charge transfer, hydrogen bonding and weak chemisorption. Therefore, the NBO 
program (Gleridening, E. D.; Reed, A E.; Carpenter, J. E., Weinhold, F. NBO Version 
3.1, 1995) was used for studying the electron density distribution of the adsorption 
system. 

Results and Discussion 

For purification, as the inventors have shown for the removal of ppm 
levels of dienes from olefins (Takahashi, A; Yang, F. H.; Yang, R. T. 
"Aromatics/Aliphatics Separation by Adsorption: New Sorbents for Selective Aromatics 
Adsorption by n -Complexation," Ind, Eng. Chem. Res., 2000, 5P, 3856\ and Takahashi, 
A; Yang, R T. "New Adsorbents for Purification: Selective Removal of Aromatics," 
AIChE J., submitted 2001), one parameter is the separation fiictor; and the separation 
factor should be evaluated at low concentrations of the solute but at high concentrations 
of the solvent In other words, the capacity for the dilute solute (at ppm level) is an 
important &ctor in determining the purification. The separation factor (a) is given by: 

ai2-XiY2/(X2Yi) (2) 

where X and Y are mole fractions in the adsorbed and fluid phases, respectively. In our 
case, component 1 is thiophene and component 2 is benzene. Hence, the ability of the 
sorbent for purification depends on the amount of thiophene adsorbed at the ppm levels. 

Figure 1 shows the equilibrium isotherms of benzene and thiophene on 
NaY at 120°C. NaY does not have a selectivity for thiophene; however, it adsorbs both 
benzene and thiophene quite strongly, as evidenced by the still measurable amounts 
adsorbed at partial pressures as low as 1 0-1 GO ppm of 1 atm. As expected, a fixed-bed 
adsorber experiment using NaY would yield a breakthrough curve (Yang, R T. Gas 
Separation by Adsorption Processes, Butterworth, Boston, 1987, Imperial College Press, 
London, 1997), as seen in the literature (Salem, A. S. H.; Hamid, H. S. "Removal of 
Sulfiir Compounds from Naphtha Solutions Using Solid Adsorbents," Chem. Eng. Tech., 
1997, 20, 342). The data were fitted by the Langmuir-Freundlich isotherm (Takahashi, 
A; Yang, RT.; Munson, C.L.; Chinn, D. "Cu(I)-Y Zeolite as a Superior Adsorbent for 
Diene/Olefin Separation," Za/igwwir, 2001; Takahashi, A; Yang, F. H.; Yang, R T. 
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"Aromatics/AIiphatics Separation by Adsorption: New Sorbents for Selective Aromatics 
Adsorption by 7C-CompIexation," Ind. Eng. Chem. Res., 2000, 39, 3856; and Takahashi, 
A; Yang, R T. "New Adsorbents for Purification: Selective Removal of Aromatics," 
AIChEJ., submitted 2001), as shown by the lines in the figure. 

Figure 2 shows the isotherms of benzene and thiophene on AgY and CuY. 
Note that: (a) compared with NaY, AgY and CuY adsorb both tiiiophene and benzene 
more strongly at below 1 0'^ atm, due to n -compiexation; and (b) on a per-cation basis, 
Cu"*" shows stronger interactions than Ag^. Compared with Figure 1, these soibents 
adsorbed significantly more thiophene and benzene than NaY at pressures below 10"^ 
atm, and nearly the same amoimts at high partial pressures. This result was a clear 
indication of tc -complexatioh with Ag**" and Cu"^; Na*^ could not form n -complexation 
bonds. From the data in Figures 1 and 2, the separation &ctors were calculated using 
mixed Langmuir-Freundlich isotherms (Taka), and are shown in Table 1 . However, these 
separation factors did not reflect the relative strengths of tt -complexation between Cu"*" 
and Ag^ because the Cu"*" exchange was not complete. Nonetheless, the large separation 
factors indicate that AgY and CuY are excellent sorbents for purification purposes. 



Table !• Separation factors of thiophene over benzene estimated from pure- 
component isotherms by extended Langmuir-Freundlich Equation at the pressure 
of 3x10'^ atm for thiophene and 1 atm for benzene. 



Adsorbent 


Ag-Y 


Cu-Y 


Na-Y 


Separation 
Factor 


2800 


50 


0.10 



The neutron activation analyses of the sorbent samples showed that the 
Ag^ exchange was 100% but the Cu"^ exchange was only 23%. On a per-cation basis, the 
n -complexation vnth Cu^ was stronger than that with Ag^. This was indeed confirmed 
by our molecular orbital calculations, as will be discussed. To understand the relative 
strengths of 7t-complexation between Ag"*^ and Cu^, the thiophene adsorption amounts at 2 
X 10'^ atm were normalized by Ag^ or Cu"*^ content, and the results are shown in Table 2. 
It is seen that Cu"^ could adsorb higher thiophene adsorption amounts per cation. In &ct, 
0.92 thiophene molecule per Cu+ was obtained at 2x10-5 atm at 120°C. This amount 
was due to Cu^ since the amount adsorbed by NaY at the same pressure was negligible. 
At the same pressure, only 0.42 thiophene/Ag**^ was obtained. This result indicated strong 
n -complexation bonds between both Cu^ and Ag"*^, and that the bond with Cu"*" was 
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Stronger. Further work is in progress to increase the Cu^ contents on various X/Y and 
other zeolites and the results will be reported shortly. 



Table 2. Thiophene adsorption amount normalized by cation content (tbiophene/cation). 



Adsoibent 


Ag-Y 


Cu-Y 


Cation for n-complexation 


Ag+ 


Cu* 


Cation Content in Zeolite (wt%) 


38.2 


8.29 


Amount Adsorbed (Qiiophaie 
molecule/cation) 


0.42 


0.92 



The heats of adsorption for u-complexation were estimated by taking the 
values at low loadings, and they were obtained from the temperature dependence of the 
isotherms at low pressures. These values are given in Table 3 . 



Table 3. Heats of adsorption Qccal/mol) for 7C-complexation. 



ZeoUte 


Thiopliene 


Benzme 


Cu(D-Y 


22.3 


21.8 


Ag-Y 


21.4 


20.1 



Bond Energies, Geometries and NBO results 

The energies of adsorption calculated using equation 1 for thiophene and 
benzene are summarized in Table 4. 



Table 4. Summaiy of enei^gies of adsoiption for tbioph^e and benzene in kcal/mol. 



MCI 


E^CTliiopliene) 


E^(Benzene) 


CuCI 


13.5 


12.4 


AgCl 


9.0 


8.6 



The theoretical calculations indicate that the n -complexation strengths follow the order 
CuCl > AgCl and more importantly, thiophene > benzene. This trend is in agreement 
with the experimental data, in Table 3. Chloride was used as the anion in the theoreticai 
calculations, while zeolite framework was the anion in the experiment It is known from 
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our previous work that the anion has a large effect on the n -complexation bonds. See^ 
for example, Huang, H,Y.; Padin, J.; Yang, R T. ''Ab Initio EfiFective Core Potential 
Study of Olefin/Paraffin Separation by Adsorption via tc -Complexation: Anion and 
Cation Eflfects on Selective Olefin Adsorption," J. Phys. Chem, B., 1999, 103, 3206; and 
Huang, HL Y.; Padin, J.; Yang, R. T. "Comparison of tt -Complexations of Ethylene and 
Carbon Monoxide and Cu"" and Ag^" Ind. Eng. Chem. Res., 1998, 38, 2720. In the 
optimized structures of thiophene-MCl complexes, the distance between the thiophene 
molecule and Cu ion is about 0.3 A shorter tiian that of thiophene and Ag ion. The NBO 
analysis is summarized in Tables S and 6. There is some donation of electron chaiges 
from the n -orbital of ihiophene to the vacant s orbital of metals known as a donation and, 
simultaneously, back donation of electron charges from the d orbitals of metals to n 
orbital of thiophene or n back-donation. The n back-donation appears to be important in 
all cases. The charge transfer results again confirmed the experimental data that the 
relative strengths of the Tc-complexation bonds follow the order thiophene > benzene and 
Cu'->Ag^ 

Table 5. Summary of NBO analysis* of n-compiexation between tbiophene and MCL 



MCI 


C-'M interaction 
(o donation) 


M-'C interaction 
(d - n* backdonation) 


Net Change 
Qi + qi 


CuCl 


0.037 


-0.022 


0.015 


AgCl 


0.022 


-0.014 


0.008 



qi is ttie amount of electron population change in valence s oibitals of the metal, and is the 
total amount of election population decrease in valence d oibitals of the metal. 



Table 6> Summary of NBO analysis* of Tt-complexation between benzene and MCL 



MCI 


C-*M interaction 
(o donation) 
Qi 


M-* C interaction 
(d - n* backdonation) 


Net Change 
qi + Oi 


CuCl 


0.011 


-0.013 


-0.002 




0.003 


-0.007 


-0.004 



* is the amount of electron population change in valence s oibitals of the metal, and is the 

total amount of dectron population decrease in valence d oibitals of £he metal. 

EXPERIMENT B 

Thiophene Removal Using Sodium and Silver Cation in Faujasite Zeolites 

Work has been done on diflferent adsorbents to remove sulfur from 
mixtures, such as, for example, thiophene from benzene or n-octane based solutions. 
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These adsorbents include sodium and silver forms of faujasite type zeolites, particularly, 
Y type zeolites. Initial studies were done for the case of an inlet concentration of 2000 
ppmw of thiophene in either benzene or n-octane. 
Sorbent Preparation 

Sodium Y Zeolite powder was obtained from Strem Chemicals. Silver Y 
Zeolite was obtained by ion exchanging NaY at room temperature. The solution was 
prepared using a silver nitrate salt and deionized water (--0.2 M). A 4-fold excess 
amount of silver was used to ensure exchange completion. After ion exchanging for 48 
hours, the adsorbent was recovered by filtration and washed v^th large amounts of 
deionized water. Drying was performed at room temperature. 
Breaktihrough Experiments 

Breakthrough experiments were performed in a custom made quartz 
reactor equipped with a glass frit The adsorbent was loaded inside the reactor (usually 1 
or 2 grams) and activated in situ under a helium atmosphere at SSO^'C for 24 hrs. The 
adsorbent was then cooled down to room temperature. At this point, the sorbent was 
wetted by pouring down sulfur free solvent (benzene or n-octane) for about 30 minutes. 
This stream was then changed to one containing 2000 ppmw thiophene. Samples were 
collected at different time intervals, and the thiophene outlet concentration was analyzed 
using a GC unit equipped with a polar column at 65**C. 
Preliminary Results 

Figure 3 shows breakthrough curves for thiophene in either benzene or n- 
octane in sodium or silver forms of Y Zeolite. 

Thiophene adsorption amounts after saturation, 
NaY AiY 
Saturation Amount (mmol/g) 

0.102 0.171 
— 0.898 



Solvent 

Benzene 
n-Octane 



Integration of the area above the curves yield the adsorption saturation 
amounts presented in the table shovm above. For the case of thiopheneA)enzene 
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mixtures, the calculated thiophene adsorption amounts are not that high. Very recent 
single component vapor phase adsorption experiments have shown that NaY and AgY 
may adsorb benzene and Ihiophene in similar quantities. Thus, benzene may compete for 
many of the adsorption sites, perhaps leaving thiophene with fewer available sites. 
However, despite the possibility of competitive adsorption, the presence of silver cations 
results in an increase in some adsorption capacity. For the case of thiophene/n-Octane 
mixtures, AgY Zeolite shows great selectivity towards thiophene. Without being bound 
to any theory, it is believed that this indicates tiiat the interaction between the silver 
cations and llie thiophene molecules is higher than that of n-Octane molecules. 
EXPERIMENT C 

In this experiment; the known commercial sorbents such as Na-Y, Na- 
ZSM5, H-USY, activated carbon and activated alumina (Alcoa Selexsorb) were included, 
and a direct comparison was made with Cu-Y and Ag-Y which were the sorbents with 7t- 
complexation capability. Thiophene and benzene were used as the model system for 
desulfurization. 
Sorbent preparation 

Various kinds of sorbents were investigated in this work. Four as-received 
sorbents: Na-type Y-zeolite (Na-Y, Si/Al=2.43, 56 Al atoms/unit cell, Strem Chemical), 
H-type ultia-stable Y-zeolite (H-USY, Si/Al=195, 0.98 AI atoms/u.c, HSZ-390HUA, 
TOSOH Corporation), activated carbon (Type PCB, Calgon Carbon Corporation) and 
modified activated alumina (Selexsorb CDX, Alcoa Lidustrial Chemical), were used in 
this study. According to the product datasheets, Selexsorb CDX is formulated for 
adsorption of sulfur-based molecules, nitrogen-based molecules, and oxygenated 
hydrocarbon molecules. Na-Y and H-USY were in powder form (binderless). Since 
activated carbon was in granular form and activated alumina was in pellet form, they 
were crushed into powder form for evaluation. 

Cu'*'-Y and Ag-Y were prepared as in Experiment A, above. 13X 
(Si/Al=1.25, Linde, lot#945084060002) was used for the preparation of Cu-X (10 fold 
CEC solution of Cu(N03)2, ion-exchanged at 65 ^'C for 24 hrs, three times) and Ag-X (5- 
fold CEC solution of AgNOa, ion-exchanged at RT for 24 hrs, twice). Na-type ZSM-5 
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(Na-ZSM-5) was prepared at room temperature by Na^-exchange of NH4-ZSM-5 
(Si/Al=10, SM-24, ALSI-PENT A Zeolite GmbH). 
Isotherm and uptake rate measurements 

The objective of this study is to compare the strength of adsorptive 
interaction between adsorbents and thiophene/benzene. Extremely low partial pressures 
at less than 10*^ atm would be necessary to meet this objective if isotherms were 
measured at ambient temperature, because the isotherms at ambient temperature are &irly 
flat and are difficult to compare each other. However, it is very difficult to obtain and 
control such low partial pressures experimentally. Therefore, single component 
isotherms for benzene and thiophene were measured at 90°C^ 120^C and 180 using 
standard gravimetric methods. A Shimadzu TGA-SO automatic recording microbalance 
was employed. Helium (Pre-purified gmde, Metro welding, 99.995%) was used as the 
carrier gas and was first passed through two consecutive gas-wash bottles (to ensure 
saturation), which contained benzene (EIPLC grade, Aldrich, 99.9+%) or thiophene 
(Aldrich, 99+%). After diluting the concentration to the desired value by blending with 
additional helium, the mixture was directed into the microbalance. 

Isosteric heats of adsorption were calculated using the Clausius-Clapeyron 
equation from isotherms at different temperatures. Nitrogen isotherms at 77K measured 
with a Micromeritics ASAP 2010 system were used for pore size distribution and pore 
volume determination. Pore size distributions were calculated with Horvath-Kawazoe 
equation. See Horvath, G., and K. Kawazoe, "Method for Calculation of Effective Pore 
Size Distribution in Molecular Sieve Carbon,*' J. Chem. Eng. Japan, 1983, 16, 470; Saito, 
A, and H. C. Foley, "Curvature and Parametric Sensitivity in Models for Adsorption in 
Micropore,'' >4/C/i£; J, 1991, 37, 429; and Cheng, L. S., and RT. Yang, "Improved 
Horvath-Kawazoe Equations Including Spherical Pore Models for Calculating Micropore 
Size Distribution," Chem. Eng. Set 1994, 49, 2599. Slit pore model was used for 
activated carbon, and cylindrical pore model was used for ZSM-S and activated alumina, 
while spherical pore model was applied to Y-zeolites and X-zeolites. As for the 
parameters for H-K equation such as polarizability, magnetic susceptibility, density, etc., 
de&ult values in Micromeritics ASAP system were used for the calculation of pore sizes. 
See "Accelerated Surface Area and Porosimetry System Operator's Manual V3.02," 



20 



wo 03/020850 



PCT/US02/28040 



Micromeritics Instrument Corp., (1997). The diffiision time constants, D/r^ ( s"^ ), were 
calculated from the uptake rates. See Yeh, Y. T, "Diffusion and Adsorption of Gases in 
Molecular Sieves," Ph.D. Dissertation, University of New York at Buffalo, Bufl&lo, New 
York, 1989. In this experiment, short time region (up to 30 % uptake) and spherical 
adsorbent model were used. 
Chemical Analysis 

The compositions of Cu-Y and Ag- Y were characterized using neutron 
activation analysis (NAA) in the research nuclear reactor of the Phoenix Memorial 
Laboratory at the University of Michigan. The sample was irradiated sequentially for one 
minute at a core-&ce location with an average thermal neutron flux of 2 x 10^^ n/cm^/s. 
Two separate gamma-ray spectra were then collected for each sample witii a high 
resolution germanium detector: one after 13-minute decay to determine the 
concentrations of Al, Ag and Cu, and a second after a 1 hour and S6-minute decay to 
analyze for Na Gamma energy lines at 1779 keV, 632.99 keV, 1039.20 keV and 1368.6 
keV were used for the determination of Al, Ag, Cu and Na concentration, respectively. 

Molecular Orbital Computational Details and Natural Bond Orbitals are as 
described in Experiment A above. 
Models for Ag-Zeolite (AgZ) and Cu-Zeolite (CuZ) 

The zeolite models selected for this study are similar to the ones used by 
Chen and Yang (see Chen, N. and RT. Yang, "Ab Liitio Molecular Orbital Study of 
Adsorption of Oxygen, Nitrogen, and Ethylene on Silver-Zeolite and Silver Halides," 
Ind, Eng Chem. Res. 1996, 35, 4020), with the molecular formula of (HO)3Si- O - 
Al(OH)3, and the cation Ag+ or Cu+ sits 2 - 3 A above the bridging oxygen between Si 
and Al. This is a good cluster model representing the chemistry of a univalent cation 
bonded on site II (SH) of the feujasite fiamewoik (Z). Once the optimized structures of 
AgZ and CuZ are obtained at the B3LYP/LanL2DZ level, then a molecule of thiophene 
(C4H4S) or benzene (CeHe) is added onto the cation of the zeolite model, and the 
resulting structure is further optimized at the B3LYP/LanL2DZ level. 
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Results and Discussion 

Characterization of sorbents. Neutron activation analysis showed that 
Ag/Al and Na/Al ratios in Ag-Y were 1.13 and 0.01, respectively. This is because Ag^ is 
known to have higher selectivity to cation sites in zeolites compared to Na"**. See Breck, 
D.W. Zeolite Molecular Sieves: Structure, Chemistry and Use\ Wiley: New York, 1974. 
More than 100% Ag ion-exchange ratio was obtained in Ag-Y, because some Ag ions 
were located outside the charge-compensating sites. On the other hand, the Al, Cu, and 
Na contents in Cu-Y were 6.10 ± 0.25 wt%, 6.65 ± 0.05 wt%, and 1.50 ± 0.12 wt%, 
respectively. This means that the Cu/Al molar ratio was 0.463, and the Na/Al ratio was 
0.289, indicating less selectivity than Ag"^. Therefore, Cu ion-exchange for X-zeolite was 
carried out three times at 65 °C (instead of twice at RT for Cu-Y) for 24 hrs to enhance 
the ion-exchange. 

Figure 4 shows the comparison of the cumulative pore volumes of the 
sorbents used. Activated carbon has the largest pore volume with small pores, v^ile Na- 
ZSM-5 has the smallest pore volume. Selexsorb CDX also shows small poies. Because 
silver is heavier than copper and the silver content in Ag-Y (Ag/Al=l . 1 3) was more than 
copper content in Cu-Y (Cu/Al = 0.463), the density of Ag-Y was higher by 36% 
0^ao.o84Cuo.i35Aloj»2Sio.70802; MW = 70.27 vs. Ago.33aAlo.292Sio.70802; MW = 95.36). 
Consequently, the pore volume of Ag-Y was about 25 % smaller than that of Cu-Y. 

Benzene/Thiophene adsorption isotherms: Na^Y. Figure 5 shows the 
equilibrium isotherms of benzene and thiophene on Na-Y at 120°C and ISO^C. The data 
were fitted by the Langmuir-Freundlich and Dubinin-Astakhov isotherms (see Padin, X; 
Yang, R. T.; Munson, C. L., "New Sorbents for Olefin-ParaflBn Separations and Olefin 
Purification for C4 Hydrocarbons," Ind. Eng. Chem. Res., 1999, 38, 3614; Jayaiaman, A; 
Yang, R T.; Munson, C. L.; Chinn, D., "Deactivation of 7c-Complexation Adsorbents by 
Hydrogen and Rejuvenation by Oxidation," Ind. Eng. Chem. Res., 2001, 40, 4370; 
Takahashi, A.; Yang, R. T.; Munson, C. L.; Chinn, D., "Cu(r)-Y Zeolite as a Superior 
Adsorbent for Diene/Olefin Separation," Z,ang7wi/z>, 2001, 17, 8405; Takahashi, A.; 
Yang, F. H; Yang, R T., "Aromatics/Aliphatics Separation by Adsorption: New 
Sorbents for Selective Aromatics Adsorption by Tr-Complexation," Ind. Eng. Chem. Res., 
2000, 39, 3856; and Takahashi, A.; Yang, R T, "New Adsorbents for Purification: 
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Selective Removal of Aromatics " AIChEJ., submitted (2001)), which are shown by the 
lines in Figure S. Barthomeuf and Ha measured benzene adsorption on NaH-Y 
(Si/Al=2.43, 81%Na) at 172 See, Barthomeuf, D. and Ha, B.-K, "Adsorption of 
Benzene and Cyclohexane on Faujasite-type Zeolites, Part-1 , Thermodynamics Properties 
at Low Coverage," J. Chem. Soc, Faraday Trans. 7, 1973, 69, 2147. Their isotherm and 
calculated heat of adsorption agreed with the data shown in Figure 5 very well, 
substantially ensuring the accuracy of the isotherms in this work. More benzene was 
adsorbed on Na-Y than thiophene at pressures lower than 1 0'^ atm. This result indicated 
that the contribution of the higher polarizability of benzene (10.3x10**^ cm^ for benzene 
vs. 9.7x10"^^ cm^ for thiophene) was larger than the contribution of the dipole moment (0 
debye in benzene vs. 0.55 debye in thiophene) and lower magnetic susceptibility (9.1x10" 
^ cm^/molecules for benzene vs. 9.5x10"^ cm^/molecules for thiophene). At higher 
pressures over 1 x 10"^ atm, the benzene adsorption amount became smaller owing to the 
larger molar volume of benzene than thiophene. Na-Y does not have a selectivity for 
thiophene; however, it adsorbs both benzene and thiophene quite strongly, as evidenced 
by the still measurable amounts adsorbed at partial pressures as low as 10"^ to 1 0'^ atm. 

Ag-Yand Cu^Y. Figure 6 and Figure 7 show the isotherms of benzene and 
thiophene on Ag-Y and Cu-Y. Compared with Figure 5, these sorbents adsorbed 
significantly more thiophene/benzene than Na-Y at pressures below 1 0'^ atm, and nearly 
the same amounts at high partial pressures. This result was a clear indication of tc- 
complexation with Ag^and Cu"*"; since Na"^ could not form 7c-complexation bonds. 
However, the difference in the amounts of thiophene/benzene adsorbed did not reflect the 
relative strengths of 7i-complexation between Cu"*" and Ag*^ because the Cu"*" exchange was 
not complete. The neutron activation analyses of the sorbent samples showed that the 
Ag"*^ exchange was 100% but the Cu"*" exchange was only 46%. According to the EPR 
analysis (described in Takahashi, A., Yang, R. T., Munson, C. L. and Chinn, D., "Cu(I)- 
Y Zeolite as a Superior Adsorbent for Diene/Olefin Separation," Langmuir, 2001, 17, 
8405), only about a half of the Cu^^ was auto-reduced to Cu^ after our heat treatment at 
450°C for Ihr in He. On a per-cation basis, the Tc-complexation with Cu"*^ was stronger 
than that with Ag^. This was indeed confirmed by our molecular orbital calculations, as 
will be discussed. To understand the relative strengths of 7i-complexation between Ag^ 
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and Cu"^, the thiophene adsorption amounts at 2 x 10"^ atm were normalized by Ag"^ or 
Cu^ content, and the results are shown in Table 7. It is seen that Cu^ could adsorb higher 
thiophene adsorption amounts per cation. In feet, 0,92 thiophene molecule per Cu^ was 
obtained at 2x10'^ atm at 120^C. This amount was due to Cu^ since the amount adsorbed 
by NaY at the same pressure was negligible. At the same pressure, only 0.42 
thiophene/Ag^ was obtained. This result indicated strong Tc-complexation bonds between 
both Cu^ and Ag^, and that the bond with Cu**^ was stronger. 

The vapor-phase isotherms of benzene/thiophene on Cu-Yand Ag-Y were 
not completely reversible within a practical time scale (1 hours to 3 hours), because of 
relatively strong interaction. However, the bond energies are of the order of 20 kcal/mol, 
which are weak enough so a substantial fiaction of the adsorbed molecules are desorbed 
in a short period of time. For practical application, thermal desorption or displacement 
desorption may be a desirable option. 

Cw-X In order to increase the Cu^ content in the feujasite structure, an X- 
zeolite (Si/Al = 1 ,25) was ion-exchanged with Cu(N03)2. Cu-X was considered to be 
proinising, since larger amounts of Cu"*" can be present in the structure of zeolite, 
especially at SIU sites. The color of Cu-X powder after ion-exchange (before auto- 
reduction) was bluer than that of Cu-Y, suggesting that more Cu^^ was contained in X- 
zeolite than Cu-Y. For faujasite-type zeolites, the cations are designated as SI (the center 
of hexagonal prism), ST (opposite SI but located in the cubooctahedron), Sn (single six- 
ring in the supercage), SII* (opposite Sn but inside the cubooctahedron), and SDI (near 
the four-ring windows in the supercage). See Yang, R. T. Gas Separation by Adsorption 
Processes, Butterworth, Boston, 1987, Imperial College Press, London, 1997; and 
Hutson, N. D.; Reisner, B. A.; Yang, R. T.; Toby, B. H., "Silver Ion-Exchanged Zeolites 
Y, X, and Low-Silica X: Observations of Thermally Induced Cation/Cluster Migration 
and the Resulting Effects on the Equilibrium Adsorption of Nitrogen," Chem. Mater, 
2000, 12, 3020. The cations at SI, SF and SIT sites are not exposed to the supercage and 
are shielded by framework oxygen. They cannot interact with molecules inside the 
supercages directly. Only SH and SIU sites are exposed to the supercage and can interact 
with molecules inside the supercage. Because Y-zeolite (Si/Al=2.43) does not have the 
Sm cation sites, it is believed that only cations at SH site are responsible for the 
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thiophene/benzene adsorption. On the other hand, X-zeolite (Si/Al=l .25) have a certain 
amount of SIDE sites, based on the results on Ag-X. This additional Cu^ located at SHI 
site should enhance the adsorption amount However, unfortunately, it was found that 
Cu-X (Si/Al=l .25) was not stable at temperatures over 200°C, and the pore volume was 
decreased by 60 %. The color of the powder was changed fix>m light blue to dark 
greenish gray after heating at 200 X for 1 hr in vacuum. As a result, thiophene 
adsorption amount on Cu-X after 450 auto-reduction was considerably smaller than 
that on Cu-Y (Isotherms were not shown). To alleviate the destruction of Cu-X structure, 
the heating rate to 450 °C was reduced to 0.5°C/min and 2''C/min (instead of 10 ^'C/min ), 
but little improvement was observed. Although the reason for poor stability is not clear 
yet, copper cations were responsible for this phenomenon, because Ag-X, which will be 
discussed next, did not show any degradation of pore structure. It was reported that 
Nickel exchanged A-zeolite was unstable and lost its crystal structure at a temperature 
above 70 °C. See Breck. D.W., Zeolite Molecular Sieves: Structure, Chemistry and Use; 
Wiley: New York, 1974. The thermal instability of nickel exchanged A-zeolite was 
interpreted in terms of a ligand-field effect After dehydration octahedral coordination 
for the nickel cation is not &vored in the zeolite structure, A similar phenomenon might 
be occurred in copper exchanged X-zeolite. 

Ag'X, Thiophene adsorption performance on Ag-X was also examined to 
imderstand the effect of additional Ag^ at SHI sites. Thiophene adsorption amoimts on 
Ag-Y and Ag-X at 120°C were compared in Figure 8. Contrary to our expectation, Ag-X 
adsorbed less thiophene than Ag-Y on a per weight basis. However, since Ag-Y and Ag- 
X have different densities due to their different compositions, the adsorption amounts of 
thiophene were nearly the same based on per unit cell of ftiujasite (i.e., about 25 
tiiiophene molecules/u.c. at 2.3x1 0"^ atm and 36 — 40 thiophene molecules/u.c at 1 .8x10"^ 
atm for both Ag-Y and Ag-X). The reason for the ineffectiveness of extra Ag^ at SID site 
in Ag-X is not known. The thiophene adsorption amounts were also examined fiom the 
viewpoint of zeolite pore volume and liquid thiophene volume. The molar density of 
liquid thiophene is calculated to be 83 cc/mol at 120 °C. See Doong, S. J.; Yang, K T., 
"A Simple Potential-Theory Model for Predicting Mixed-Gas Adsorption,*' Ind. Eng. 
Chem. Res. 27, 630. 1.5 mmol thiophene adsorbed per gram of Ag-Y at 2.3x10'^ 
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atm corresponded to 0.125 cc liquid thiophene per gram of Ag-Y. The thiophene liquid 
volume of 0. 125 cc was nearly one half of pore volume of Ag-Y (0.24cc/g as shown in 
Figure 4), suggesting that one half of the pore was filled with thiophene. Therefore, it 
does seem likely that the cavity of Ag-X was filled with thiophene even at 2.3x10"^ atm 
and only a small amount could be further adsorbed at higher pressures. One possible 
reason for the relatively low amounts in Ag-X is that thiophene molecules adsorbed could 
be located at positions near more than one silver cations at both SR and SM sites 
simultaneously, so that no improvement in adsorbed amount was observed. 

Silver in Ag-X is known to form silver cluster at elevated temperature. 
See Sun, T. and K Sefl^ "Silver Clusters and Chemistry in Zeolites," Chem. Rev. 1994, 
94, 857. Hutson et al. showed that nitrogen adsorption capacity on Ag ion-exchanged 
low-silica X-zeolite (Si/Al =1) was enhanced by heating in vacuum at 450°C for more 
than 4 hours. This enhancement was explained by thermally induced Ag cation and/or 
cluster migration using Rietvelt refinement of neutron powder diffraction data. The 
thiophene adsorption isotherm on Ag-X(Si/Al = 1 .25) after 450 ^C for 4 hours is also 
shown in Figure 5. No enhancement was observed for thiophene adsorption on Ag- 
X(Si/Al=1.25). 

H'USY. Thiophene and benzene isotherms on H-USY(Si/Al = 195) are 
shown in Figure 9. The interaction between thiophene/benzene and high-silica H-USY 
was very weak. Neither thiophene nor benzene was adsorbed at pressures lower than 
1x10"^ atm. Adsorbed amounts were substantially smaller than Na-Y even at 2x1 0"^ atm. 
These results are clearly due to the strong interactions between the cations in zeolite and 
the thiophene or benzene molecule, 

Na-ZSM-S. As stated above, several groups used Na-ZSM-5 for 
thiophene/benzene purification and showed that Na-ZSM-5 could remove thiophene 
impurities fi*om benzene in their fixed bed breakthrough experiments. Figure 10 shows 
isotherms on Na-ZSM-5. Although thiophene and benzene isotherms turned out to be 
virtually the same on Na-ZSM-5, it may be important to note that thiophene isotherms 
were feirly flat over the pressure range in the application of purification. This small 
difference of adsorbed amounts at 3x10"^ atm and 2x10'^ atm was preferable for the 
purification of benzene by removal of thiophene, which will be discussed below. One 
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disadvantage of Na-ZSM-5 was small adsorption capacity owing to the small pore 
volume, as clearly seen in Figure 4. 

Activated carbon and modified activated alumina. The thiophene and 
benzene adsorption isotherms on activated carbon (Type PCB) and modified activated 
alumina (Selexsorb CDX) are shown in Figure 1 1 and Figure 12. PCB activated carbon 
is commercially designated for use in liquid phase or vapor phase applications such as 
recovery of alcohols, chlorinated hydrocarbons, hydrocarbons and aromatics. Selexsorb 
CDX is specially fonhulated by Alcoa Industrial Chemicals for adsorption of polar 
organic compounds including sulfur-based molecules (mercaptans, sulfides, disulfides, 
thiophenes), nitiogen-based molecules (nitriles, amines, pyridines) and oxygenated 
hydrocarbon molecules (alcohol, glycols, aldehydes, ketons, ethers, peroxides). More 
benzene than thiophene was adsorbed on both activated carbon and modified activated 
alumina at pressures below 1x10'^ atm. This is surprising because Selexsorb was thought 
to adsorb thiophene more selectively than benzene. This result indicates that 
polarizability may play a more important role (than permanent dipole moment and 
magnetic susceptibility) in adsorption on these two sorbents. At pressures higher than 1 x 
10"^ atm, pore filling dominated hence benzene adsorbed less than thiophene because of 
the larger molar volume of benzene than thiophene. 

Comparison of thiophene adsorption on all sorbents, Thiophene 
adsorption isotherms on all sorbents are compared in Figure 13. It is clearly seen that 
Ag-Y and Cu-Y could adsorb significantly larger amounts of thiophene even at very low 
pressures. 

Purification capabUity 

For purification, as has been shown for the removal of ppm levels of 
dienes firom olefins, an important consideration is the separation factor; and the 
separation fiictor should be evaluated at low concentrations of the solute but at high 
concentrations of the solvent In other words, the capacity for the dilute solute (at ppm 
level) is a &ctor that determines the purification. As stated above, the separation factor 
(a) is given by: 

ai2 = XiY2/(X2Y0 (2) 
where X and Y are mole fractions in the adsorbed and fluid phases, respectively. In our 
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case, component 1 is thiophene and component 2 is benzene. A further important 
parameter is the capacity of the sorbent for the solute at low concentrations. Hence, the 
ability of the sorbent for purification depends on the amount of thiophene adsorbed at the 
ppm levels. 

To estimate the adsorption amount in the mixture of thiophene and 
benzene, the pure-component adsorption data were first fitted by both D-A (Dubinin- 
Astakhov) and L-F (Langmuir-Freundlich) isotherms. The separation fiictors of 
thiophene over benzene were then obtained using the equivalent multi-component 
isotherms, as described previously. The conversion of molar adsorbed amount to 
volumetric adsorbed amount and vice versa were performed in the same way as Doong 
and Yang. The molar volumes of thiophene in this work were 79.7 cc/mol at 90 °C, 83.0 
cc/mol at 120*^0 and 89.8 cc/mol at 180°C, and those of benzene were 90.5 cc/mol at 
90°C, 94.6 cc/mol at 120°C and 102.9 cc/mol at 180*^0. The adsorbed molar volumes 
were assumed to be the same in both single- and mixed-gas adsorption. Separation 
fiictors were calculated for two cases: lxlO"^atm thiophene/1 atm benzene and 1x10"^ atm 
thiophene/1 atm benzene. 

The parameters for the D-A and L-F equations are given in Table 8 and 
the resulting separation factors are listed in Table 9. Some of the D-A parameters were 
not applicable for the Doong/Yang model for mixtures, when the volumetric adsorbed 
amount of benzene was much higher than that of thiophene. The separation fectors of 
thiophene over benzene follow the order: Ag-Y > Na-ZSM-5 > Cu-Y = activated carbon 
» Na-Y » H-US Y = modified activated alumina. This order can be explained as 
follows. As described earlier, the ability of the sorbent for purification depends on the 
amount of thiophene adsorbed at the ppm levels, compared with the amount of benzene 
adsorbed at high concentration. Ag-Y and Na-ZSM-5 showed feirly flat isotherms of 
thiophene over the entire partial pressure range, resulting in the higher separation &ctors. 
In feet, Ag-Y and Na-ZSM-5 maintained 65% and 47% of thiophene adsorption amounts 
at 2.3x10"^ atm compared with those at 1 .3x10"^ atm, respectively. The advantage of Ag- 
Y against Na-ZSM-5 was fiirther enhanced by the fact that Ag-Y was slightly selective 
for thiophene than benzene, while thiophene/benzene isotherms on Na-ZSM-5 were 
virtually the same. On the other hand, Cu-Y and activated carbon could adsorb a certain 



28 



wo 03/020850 



PCT/US02/28040 



amount of thiophene/benzene, but smaller adsorption amounts at lower pressures reduced 
the separation Actors considerably. Only 1 8% and 5% of thiophene adsorption amounts 
were retained by Cu-Y and activated carbon between 2.5 x 1 0" atm and 1.3x10* atm, 
respectively. Larger adsorption amounts of benzene than thiophene by Selexsorb CDX 
made the separation &ctor of this sorbent extremely low. 

Another important parameter for purification is the sorbent capacity for 
the solute at very low concentrations, i.e., for thiophene at ppm levels. The sorbent 
capacities are shown in Figure 13. At 2.3x10'^ atm, the capacities for thiophene follow 
the order Ag-Y & Cu-Y » Na-ZSM-5 > activated carbon > Na-Y > modified alumina 
& H-USY. Based on these two key parameters, the Tc-complexation sorbents are clearly 
the best 

Heat of adsorption 

Heats of adsorption were calculated using the Clausius-Clapeyron 
equation fi-om isotherms at two different temperatures, and are shown in Table 10. All 
the heats of thiophene/benzene adsorption had the tendency to decease as the loading 
increased. This is a common phenomenon for the sorbents such as ion-exchanged 
zeolites that have heterogeneous sites. The heats of adsorption on activated carbon, in 
particular, ranged widely fi-om 23.9 kcal/mol (at O.Smmol/g loading) to 8.0 kcal/mol (at 
3mmol/g). Ag-Y and Cu-Y exhibited the higher heats of adsorption than Na-Y for both 
benzene and thiophene; it is believed due to Tc-complexation. More importantly, the heats 
of adsorption for Aiophene were higher than that of benzene. These experimental results 
can be explained by molecular orbital calculation and NBO analysis, which will be 
discussed shortly. 

At the low loading of O.S mmol/g, Na-ZSM-S showed nearly the same 
heats of adsorption as Na-Y for both thiophene and benzene. The different pore 
dimensions (5.2-5.6 A for ZSM-5 vs. 7.4 A for Na-Y) apparently had no influence on the 
heats of adsorption. It is not clear why the amounts adsorbed on Na-Y decreased sharply 
at very low pressures while that on Na-ZSM-5 maintained. 
Bond Energies, Geometries and NBO results 

The energies of adsorption calculated using equation 1 for thiophene and 
benzene are summarized in Table 1 1. The theoretical calculations indicate that the 
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complexation strengths follow the order CuZ > AgZ and more importantly, thiophene > 
benzene. This trend is in agreement with the experimental data, in Table 10. In fact, the 
molecular orbital results on CuZ and AgZ are in excellent agreement with the 
experimental data. Both chloride and zeolite models were used as the anion in the 
theoretical calculations, while only zeolite fiamework was the anion in the experiment It 
is believed that the anion has a large effect on the 7C-complexation bonds. The bond 
energies on the zeolites (Z) are significantly higher than that on the chlorides (Table 1 1). 
This result indicates that the zeolite anion is more electronegative than the chloride anion, 
which has been revealed by Chen and Yang in their ab initio molecular orbital 
calculations. In the optimized structures of thiophene-MCl complexes, the distance 
between the thiophene molecule and Cu ion is about 0.3 A shorter than that of thiophene 
and Ag ion for chloride. The NBO analysis is summarized in Tables 12 and 13. There is 
some donation of electron charges from the 7C-orbital of thiophene to the vacant s orbital 
of metals known as a donation and, simultaneously, back donation of electron charges 
from the d orbitals of metals to n orbital (i.e., anti-bonding n orbital) of thiophene om 
back-donation. It appears that the a donation is more predominant for thiophene and the 
n back-donation is more important for benzene. Comparing the two anions, zeolite anion 
and chloride anion, the NBO results show that both a donation and d - tc* backdonation 
are significantly stronger with the zeolite anion bonded to Ag^ or Cu"*^. The charge 
transfer results again confirmed the experimental data that the relative strengths of the tc- 
complexation bonds follow the order: thiophene > benzene and Cu^ > Ag^ 
Uptake Rates 

Table 14 compares difiiision time constants for Ag-Y, Cu-Y and Na-Y. 
These constants were calculated from the uptake rates using the solution for Pick's law. 
Uptake rates among the three soibents were very similar and no notable differences 
between thiophene and benzene were observed. These rates are reasonably &st for 
practical applications. 

In this experiment, vapor-phase benzene/thiophene adsorption isotherms 
were investigated to develop new sorbents for desulfurization. Among the sorbents 
studied, it was found that Cu-Y and Ag-Y exhibited excellent adsorption performance 
(capacities and separation fectors) for desulfiirization. This enhanced performance 
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compared to Na-Y is believed to be due to the Ti-complexation with Cu"^ and Ag"**. 
Molecular orbital calculations confirmed the relative strengths of 7t-complexation: 
thiophene > benzene and Cu^ > Ag^ In an actual desulfurization process, liquid phase 
adsorptions such as batch adsorption or fixed-bed adsorption column around room 
temperatures are viable options. Liquid phase experiments are described herein. Even in 
the case of liquid phase separation, enhanced selectivities based on 7c-complexation 
studied in this experiment are applicable. 

TableT. Thiophene Adsorption Amount Normalized by Cation Content, 
(thiophene/cation molar ratio) 



Adsorbent 


Ag-Y 


Cu-Y 


Cations for n-complexation 


Ag* 


Cu* 


Cation Content in Zeolite (wt%) 


38.2 


8.29 


Amount Adsorbed (molar ratio) 


0.42 


0.92 
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3 Parameters for the P-A and L-F Equations. 



L-F Equation 









qo 


B 


n 


qo 


Vo 


C 


n 


Ps 








fmmol/er'i fl/atm'i 


(A 


fmmol/g) 




(•^ 


n 


fatm') 


Ag.Y 


Benzene 


120 °C 


3.133 


3.40 


7.67 


2,416 


0.2286 


0.0559 


1.56 


2.96 






180 


2.032 


155.5 


1.85 


2.008 


0.2066 


0.0786 


4.57 


10.07 




Thiophen 


120 •C 


3.200 


4.416 


6.59 


2.986 


0.248 


0.0609 


1.12 


2.68 






180 «C 


2.602 


6.13 


AJ7 


2.329 


0.2091 


0.0699 


1.87 


9.42 


CU-Y 


Benzene 


90 


2.989 


292.6 


1.76 


2.940 


0.266 


0.0889 


4.36 


1.34 






120 "C 


2.897 


178.6 


1.59 


2.929 


0.2771 


0.1008 


3.95 


2.96 




Tliiophen 


90 «C 


3.910 


196.4 


1.80 


3.826 


0.3048 


0.0935 


3.99 


1.19 






120 ''C 


3.548 


121.5 


1.65 


3.548 


0.2946 


0.1010 


3.81 


2.68 


Na-Y 


Benzene 


120 "C 


2.781 


13985 


0.89 


2.789 


0.2639 


0.0980 


6.80 


2.96 






180 "C 


2.820 


449.0 


0.90 


2.872 


0.2956 


0.1200 


4.73 


10.07 




Thiophen 


120 


2.001 


lo54o 


0.84 


2.658 


0.2207 


0.1050 


7.30 


2.68 






180 -C 


2.739 


15001 


0.50 


2.846 


0.2555 


0.1330 


7.54 


9.42 


PCB-AC 


Benzene 


90 «C 


/.o2o 


4.47 


2.46 


6.329 


0.5725 


0.1638 


1.71 


1.34 






120 »C 


6.312 


4.254 


2.16 


6.197 


0.5863 


0.1847 


1.62 


2.96 




Thiophen 


90 °C 


/.2ol 


5.54 


2.29 


8.936 


0.7118 


0.2246 


1.24 


1.19 






120 "C 


6.996 


6.39 


1.66 


7,236 


0.6009 


0.2193 


1.59 


2.68 


NaZSMS 


Benzene 


90 


0.928 


17.15 


3.43 


0.987 


0.0893 


0.0810 


1.44 


1.34 






120 


0.795 


14.66 


3.47 


0.826 


0.0781 


0.0793 


1,64 


2.96 




Thiophen 


90 ^'C 


1.023 


16,91 


3.27 


1.036 


0.0825 


0.0880 


1,79 


1.19 






120 *»C 


0.811 


16.36 


3.30 


0.856 


0.0711 


0.0832 


1.70 


2.68 


Seiexsorb CDX 


Benzene 


90 ^'C 


1.182 


1196,6 


1.39 


1.176 


0.1064 


0.0914 


5.23 


1.34 






120 °C 


1.096 


75902 


0.71 


1.109 


0.1049 


0.1049 


8.14 


2.96 




Thiophen 


90 ''C 


1.223 


29720 


0.79 


1.239 


0.0987 


0.1141 


6.78 


1.19 






120 "C 


1.134 


133016 


0.55 


1.158 


0.0962 


0.1265 


8.20 


2.68 


H-USY 


Benzene 


90 


4.258 


568.1 


0.48 


4.359 


0.3943 


0,2800 


4,21 


1.34 






120 "C 


3.895 


63.573 


0.53 


4.100 


0.3879 


0.2825 


3.63 


2.96 




Thiophen 


90 


4.568 


382.1 


0.51 


4.689 


0.3735 


0.2939 


3.55 


1.19 






120 "C 


4.356 


24.41 


0.61 


4.599 


0.3819 


0.2820 


4.07 


2.68 



D-A Equation 
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Table Separation Factors estimated by the D-A and L-F Equations at 120 ^C. 





LF Equation 
lE.4/1 lE-5/1 


DA Equation 
lE-4/1 lE-5/1 


AgY 


3300 


23000 


1600 


11000 


Cu-Y 


31 


78 


360 


700 


Na-Y 


0.16 


0.10 


• «»# 1) 


*** 1) 


PCBAC 


64 


160 


33 


83 


NaZSMS 


700 


3500 


1) 


1) 


Selexsorb CDX 


9.6E-04 


1.5E-04 


««* 1) 


•*» 1) 


H-USY 


1.2E-03 


2.7E04 


*•» 1) 


*** 1) . , 



1) Not applicable to Doong/Yang extension. 
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Table > Heat of adsorption (kcal/moD calculated from isotherms at different temperatures. 





Na-Y Ag-Y Cu-Y H-USY Na-2SM5 


Activated Carbon 


Activated Alumina 




(Si/AI=2.43) {SiA«M=2.43) (Si/AI=2.43) (Si/AI«195) (Si/AI=10) 


Type PCB 


Selexsorb CDX 


Thiophene 


19.1-19.6 21.3-21.5 20.8-22.4 7.9-11.2 18.6-19.2 


8.0-23.9 


16.1-17.5 




(0.5-2.0) (1.5-1.7) (2.0-3.0) (0.1-0.3) (0.45-0.60) 


(0,5-3.0) 


(0.2-1.0) 


Benzene 


17.0-18.2 19.0-20.1 19.3-21.8 6.6-13.1 16.5-17.9 


13.1-16.1 


16.8-19.6 




(1.5-2.0) (1.5-1.8) (1.8-2.5) (0.1-0.5) (0.45-0.65) 


(1.0-3.0) 


(0.6-1.0) 



*) Numbers in parentheses indicate the adsorption amounts (mmol/g) for calculation. 



The heats of adsorption decreased with loading in all cases. 
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Table 1 1 . Summary of energies of adsorption for thiophene and benzene in kcal/mol 
calculated from molecular orbital theory (Z denotes Zeolite anion) 



Eads(Thiophene) Eads(Benzene) 



CuCl 13.5 12.4 

AgCl 9.0 8.6 

CuZ 21.4 20.5 

AgZ 20.0 19.1 



Table 12. Summary of NBO analysis* of Tc-complexation between thiophene and 
MCl/MZ 

C-^M interaction M->C interaction 

(a donation) (d - n backdonation) Net Change 



MX qi q2 qi-^"q2 

CuCl 0.037 -0.022 0.015 

AgCl 0.022 -0.014 0.008 

CuZ 0.112 -0.063 0.049 

AgZ 0.101 -0.086 0.015 



* qi is the amount of electron population change in valence s orbitals of the metal, and qa 
is tiie total amount of electron population decrease in valence d orbitals of the metal. 



Table 13. Summary of NBO analysis* of 7i-complexation between benzene and MCl/MZ 

C— >M interaction M->C interaction 

(a donation) (d - 7c* backdonation) Net Change 

MX qi q2 qi+q2 

CuCl 0.011 -0.013 -0.002 

AgCl 0.003 -0.007 -0.004 

CuX 0.078 -0.119 -0.041 

AgX 0.043 -0.045 -0.002 

* qi is the amount of electron population change in valence s orbitals of the metal, and q2 
is the total amount of electron population decrease in valence d orbitals of the metal. 
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Table IH Diffusion Time Constants (1/s) of Thiophene and Benzene 





ThioDhene 




Benzene 






Pressure Change 


D/r^ 


Pressure Change 


D/r2 




(atm) 


Cl/s) 


(atm) 


(1/S-) 


Ag-Y 


1.9E-5 - 1.6E-4 


2.0E05 


5.7E-5 ■ 5.5E-4 


5.4E-05 




5.6E-4 - 1.8E-3 


3.3E-04 


5.5E-4 - 5.4E-3 


1.8E-04 




1.8E-3 ■ 1.5E-2 


1.2E-04 


5.4E-3 - 2.6E-2 


4.2E-04 


Cu-Y 


2.3E-5 - 8.8E-5 


1.3E-05 


3.2E-5 - 1.2E-4 


2.0E-05 




8.8E-5 - 4.0E-4 


2.3E-05 


1.2E-4 - 5.4E-4 


5.2E-05 




4.0E-4 - 2.2E-3 


1.5E-04 


5.4E-4 - 3.0E-3 


3.3E-04 




2.2E-3 ■ 1.3E-2 


1.4E-04 


3.0E-3 ■ 1.8E-2 


3.7E-04 


Na-Y 


8.8E-5 - 4.0E-4 


1.7E-05 


0 - 1.5E-4 


1.5E-05 




4.0E-4 ■ 2.2E-3 


1.6E-04 


1.5E-4 ■ 5.5E-4 


2.6E-05 




2.2E-3 . 1.3E-2 


2.8E-04 


5.5E-4 • 5.4E-3 


5.6E-04 




1.3E-2- 2.9E.2 


5.4E-04 


5.4E-3 - 2.6E-2 


2.4E-04 
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As mentioned hereinabove, it has further been fortuitously found that 
the sorbents of the present invention also are excellent sorbents for removal of 
aiS)matics. A detailed description of removal of aromatics follows. 
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Adsorption of benzene and cyclohexane on various Y-zeolites was investigated in 
order to develop sorbents for the purification of aliphatics by removal of aromatics. Ag-Y showed 
superior bMzene/cyclohexane selectivities than Na-Y, Pd-Y and H-US Y. Separation &Gtors 
greater than 1 0^ were obtained vsHith Ag-Y at low concentrations of benzene. The high 
selectivities were achieved by the strong interaction between benzene and Ag-Y, vsdiile the 
interaction with cyclohexane was not influenced by the cations. Molecular orbital calculation 
revealed that benzene formed classical 7c-con4)lexation bond with Ag-Y: donation of electron 
charges from the p-orbitals of benzatie to the vacant s-orbital of the silver (a donation) and, 
simultaneously, back donation of electron charges from d-orbitals of silver to tc* orbital of 
benzene (d-7r* back donation). Grand canonical Monte Carlo simulations were also performed for 
the adsorption isothemis. Potential parameters of benzene on Ag-Y including 7C-conq)lexation 
were first developed. GCMC simulation and e?q)eriments were in excellent agreement for 
adsorption of benzene on Y-zeolites. 

Purification of aliphatics fay removing aromatics is important in both petrochemical 
industry and for pollution control Because benzene (and aromatics) is known to be highfy 
carcinogenic, concentration of benz^e needs to be minimized in the automotive fuel. A number 
of separation processes have been employed to reduce the benzene concentration in a refin^'s 
gasoline pool, so that it meets new reformulated gasoline requirements. In a typical benzene 
reduction process, a combination of extraction and distillation is used (Jeanneret, 1997). Also, 
removal of aromatics from kerosene improves the clean burning properties of the fiiel and 
separation of the aromatics from the isoparaffins and naphthenes in lubricating oils improves the 
viscosity-temperature relationship (Bailes, 1983). 

Because of the importance of aromatics/aliphatics separation and the problems 
associated with current separation processes, possible alt^natives have been under continuing 
investigatioa These include pervaporation (Hao et al, 1997), liquid membrane (Li, 1968), and 
the use of liquid inclusion complexes (Atwood, 1984). Purification of dilute aromatics from 
aUphatics (e.g., toluene and/or ^Q^lene in heptane) by temperature swing adsorption (TSA) was 
studied in the liquid phase (Matz and Knaebel, 1990). Commercially available sorbents were used: 
silica gel, activated alumina, activated carbon, zeolite 13X, and polymeric resin (XAD-7). Among 
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these sorbents, silica gel was considered the best due to its superior thermal-exchange capacity. 
However, the selectivities were low. 

Adsorption is playing an increasingly important role in separations. However, its 
utility is limited by the availability of selective sorbents. The conventional sorbents and separations 
using them are based on van der Waals and electrostatic interactions between the sorbate and the 
sorbent Studies of new sorbents using weak chemical bonds such as chemical complexation have 
begun only recently in our laboratory. As suggested King (King, 1987), chemical • 
comple)cation bonds are generally stronger than van der Waals interactions, yet weak enough to 
be reversible. Therefore, tremendous opportunities exist for developing new sorbents and new 
applications in separations by using weak chemical bonds, including various forms of 
complexation bonds. 

The TT-complexation is a sub-class of chemical complexation. It pertains to the 
main group (or d-block) transition metals, i.e., from Sc to Cu, Y to Ag, and La to Au in the 
periodic table (Cotton and Wilkinson, 1966). These metals or their ions can form the normal a- 
bond to carbon and, in addition, the unique characteristics of the d orbitals in these metals or ions 
can form bonds with unsaturated hydrocarbons in a nonclassical manner. This type of bonding is 
broadly referred to as Tc-complexation. This 7c-complexation has been seriously considered for 
olefin/paraffin separations using liquid solutions containing silver or cuprous ions (Quiim, 1971; 
Ho et al, 1988; Keller et aL,1992; Blytas, 1992; Eldridge, 1993; Safarik and Eldridge, 1988). 
These involved gas-liquid operations. While gas-solid operation can be simpler as well as more 
efificient, particularly by pressure swing adsorption, the list of attempts for developing solid tc- 
cpmplexation sorbent is a short one. 

More recently, several new sorbents based on Tc-complexation were prepared for 
selective olefin adsorption. These included Ag^-exchanged resins (Yang and Kikkinides, 1995; 
Wii et al., 1997), monolayer CuCl on pillared clays (Cheng and Yang, 1995) and monolayer 
AgNOs/SiC^ (Padin and Yang, 1997; Rege et al., 1998; Padin and Yang, 2000). Moreover, Ag 
ion-exchanged Y-Zeolite was developed to purify butene by removing trace amounts of butadiene 
(Padin et aL, 1999). The purification performance of this sorbent was superior to that of NaY 
owing to 7i-complexatioa 
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It is possible, in principle, to purify aromatics from aliphatics based on n- 
complexatioa In the benzene molecules, the carbon atom is sp^ hybridized. Hence, each carbon 
has three sp^ orbitals and another Pz orbital. The six Pz orbitals in the benzene ring form the 
conjugative % bonds. The Pz orbitals also form the antibonding orbitals, which are not 
occupied. When benzene interacts with transition metals, the 7T-orbitals of benzene could ova'lap 
with the en^ty or unfilled outer-shell s orbital of the transition metal to form a a-bond 
Furthermore, the vacant antibonding 7c*-orbitals of benzene could overlap with the d-orbitals in 
the transition metals, similar to that formed m the olefin-Cu* bond (Hiiang et al, 1999a). In this 
work, Ag^ and Pd^^ ion-exchanged zeoUtes were selected for the study of aromatics/aliphatics 
separations, since these two cations have the most promising 7c-complexation capabilities, and 
they are stable. In order to understand the effects of Tt-complexation, Ag^ and Pd^^ ion-exchajnged 
zeoUtes were compared with that of Na^ ion-exchanged zeoUte and a high-Si02 zeolite. 

Benzene and cyclohexane are an ideal pair of model compounds for studying 
selective sorbents for purificatioa These molecules have similar shapes, polarizabilities and 
volatilities (the boiling points are 80 ®C for benzene and 81 for cyclohexane). The kinetic 
diameter of benz^e, which is calculated from the minimum equilibrium cross-sectional diameter, 
is estimated to be 5.85 angstrom compared with 6.0 angstrom for cyclohexane (Breck, 1974). 
Therefore, benzene and cyclohexane were used in this work. 

In addition to the experimental investigation, molecular orbital calculations were 
performed to obtain a basic understanding of the origin of the strong interactions of Ag ion- 
exchanged zeolite with benzene. Moreover, grand canonical Monte Carlo (GCMC) simulations 
were used to simulate/predict the isotherms. Many GCMC studies have been performed for 
predicting adsorption isotherms in zeoUtes. These include N2, O2, Ar, Xe, CH4, C2H4 and C4H10 
on zeolite-A, X, Y and MFI-type Zeolite (Woods and Rowlinson, 1989; Razmus and Hall, 1991; 
Kravias and Myers, 1991; Watanabe et al, 1995). Adsorption of benzene on NaY, Heulandite and 
MFI-fype zeohtes has also been studied (Snurr et al., 1994; Kitagawa et al., 1996; Klemm et al., 
1998; Laboy et at, 1999). However, these studies did not involve the efifects of 7C-complexation 
on adsorption. In this study, benz^e adsorption isotherms on Ag ion-exchanged Y-zeolite were 
simulated 
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Experimental Details 
Sorbent Preparation 

The sorbents in this work were ion-exchanged zeolites. From earlier resiilts on 
olefin/paraflBn separations (Yang and Kikkinides, 1995; Wu et al, 1997; Cheng and Yang, 1995; 
Padin and Yang, 1997; Rege et al, 1998; Padin and Yang, 2000), Ag^ and Pd^"^ were the most 
promising ions because of thdr strong 7c-complexation abilities. 

Ion-exchange of zeolites was accomplished with aqueous solutions. Faujasite was 
selected as the zeolite framework structure for ion-exchange, because of its large pore aperture 
(7.4 angstrom), large pore volimie (0.489 ml/ml) and wide range of SiOi/AfeOs ratios(l-oo). Two 
as-received zeolites, Na-type Y-Zeolite (Na-Y, Si/Al=2.43, Strem Chemicial) and H-type ultra- 
stable Y-Zeolite (H-USY, Si/Al=195, TOSOH CORPORATION, HSZ-390HUA), were used in 
this study. Pd^^ type Y-Zeolite was prepared by exchanging Na-Y zeoUte with Pd(N03)2 
(Aldrich) in an aqueous solution (Li et aL, 1994). First, Na-Y was exchanged using 2-fold excess 
of PdCNr03)2 (0.005M) at 80 ^C for 24h. After the exchange, tiie zeolite suspaision was filtered 
and washed with copious amount of deionized water. The product was dried at 110 "^C overnight 
in air. Ag^ ion-exchange was prepared at room temperature for 24 h in the same manner as Pd^^ 
exchange, using 2-fold ^cess AgNOa (O.IM). 
Adsorption Isotherms and Uptake Rates 

Sii^le component isotherms for benzene and cyclohexane were measured using 
standard gravimetric methods. A Shimadzu TGA-50 automatic recording microbalance was 
employed. Helium (Pre-purified grade, Metro Welding, 99.995%) weis used as the carrier gas and 
was first passed through two consecutive gas-wash bottles (to ensure saturation) which contained 
benzene (HPLC grade, Sigma- Aldrich, 99.9+%) or cyclohexane (HPLC grade, Aldrich, 99.9-P/o). 
After diluting the concentration to the desired value by blending with additional helium (total 
flowrate: 250 cc/min), the mixture was directed into the microbalance. 

Isosteric heats of adsorption were calculated using the Clausius-Clapeyron 
equation fi-om isotherms at different temperatures: Nitrogen isotherms at 77K measured with a 
MiCTomeritics ASAP 2010 system were used for pore size distribution and pore volume 
detemunatioa Pore size distributions were calculated with Horvath-Kawazoe equation (Horvath 
and Kawazoe, 1983; Cheng and Yang, 1994). Pore volumes were calculated at 0.95 P/Po. The 
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overall difiiision time constants, D/r^ were calculated from the uptake rates (Yeh, 1989). In this 
work, the short time region was used. 
Molecular Orbital Calculation Details 

The Tt-complexation bonding between olefins and sorbents have been extensively 
investigated using molecular orbital (MO) calculations (Chen and Yang, 1996; Huang et aL, 
1999a, 1999b). More recCTtly, MO studies were extended to benzene adsorption on various salts 
such as CuCl and AgCl (Takahashi et al.). In this work, similar MO studies for b^izene 
adsorption on Ag ion-exchanged zeolite were performed using a cluster model of Ag-zeolite. The 
Gaussian 94 Progam (Frisch et al, 1995) in Cerius2 molecular modeling software (Bowie, 1995) 
from Molecular Simulation, Inc. was used for calculations. MO calculations for benzene and Ag 
zeolite were performed at the Hartree-Fock (HDF) and density ftmctional theory (DFT) level The 
3-21G basis set is the split-valence basis set, which has been used successfully for many 
simulations for benzene adsorption on zeolites (O'Malley and Braithwaite, 1995), protonation 
reaction of propylene and isobutene on zeolite (Viruera-Martin et al, 1993) and O2, N2 and C2HU 
adsorption on Ag-Zeolites (Chen and Yang, 1996). The reliability of this basis set has been 
con&med by the accuracy of calculation results in comparison with experimental data. Therefore, 
the 3-21 G basis set was employed for both geometry optimization and natural bond orbital 
(NBO) analysis. 

Geometry OpHmization and Energy of Adsorption Calculations 

The restricted Hartree-Fock (RHF) theory at the 3-21G level basis set was used to 
determine the geometries and the bonding energies of benzene on Ag zeolite cluster, since Ag^ 
has filled d-orbitals with spin — 1. 

The computational model used for Ag zeoUte were the same as that used by Chen 
and Yang (Chen and Yang, 1996), since this model yielded heats of adsorption of N2 and C2H4 on 
Ag zeolite that agreed well with experimental data. The simplest cluster model, (H0)3Si-0(Ag)- 
A1(0H)3, was chosen for this work. The optimized structures were then used for bond energy 
calculations according to the following expression: 

Eads ~ Eadsoibate + Eadsoibent " Eadsoibent-^idsoibate ( 1 ) 

Where Eadaoibate is the total energy of brazene, EadsoAcnt is the total energy of the bare adsorbent i.e., 
the Ag zeohte cluster, and EadsoAcnt-adsoibate is the total energy of the adsorbate/adsorbent system. 
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Natural Bond Orbital (NBO) Analysis 

The optimized structxires were also used for NBO analysis at the B3LYP/3-21G 
level. The B3LYP (Becke^ 1993a) approach is one of the most useful self-consistent hybrid 
(SCH) approaches (Becke, 1993b), which is Becke's 3-parameter nonlocal exchange functional 
(Becke, 1992) with the nonlocal correlation fimctional of Lee, Yang and Parr (Lee et al., 1988). 

The NBO analysis performs population analysis pertaining to the localized wave- 
function properties. It gives a better description of the electron distribution in compounds of high . 
ionic character, such as those containing metal atoms (Eleed et aL, 1985). It is known to be 
sensitive for calculating localized weak interactions, such as charge transfer, hydrogen bonding 
and weak chemisorption. Therefore, the NBO program (Glendening et al., 1995) was used for 
studying the electron density distribution of the adsorption system. 
Monte Cario Simulation Details 

Grand canonical ensemble (at a fixed pressure) Monte Carlo simulation was 
performed utilizing the Cerius2 molecular modeling software (Bowie, 1995) on a Silicon Graphics 
Indigo2 workstation running IRIX v. 6. 5. After the "crystal builder'' module was used to set up 
the zeolite models, energy expressions and parameters for force field were input by the "force 
field editor" module. The "sorption" module was then used for benzene adsorption simulations. 

The simulation was performed at 120 °C and 180 for at least 1,000,000 
configurations. During the simulation, the sorbate molecules interacted with the potential field 
generated by sorbate-sorbent and sorbate-sorbate interactions. The initial configurations contained 
no sorbate molecules. Each subsequent configuration was generated by one of four moves: 
CREATE, DESTROY, TRANSLATE and ROTATE. Molecular creation attempts were made at 
random points within the accessible portion of the zeolite lattice with the criterion that there be no 
overkpping sites (i.e. creations which results in interaction sites which were closer than half the 
sum of the van der Waal's radii of the two sites were rejected). Translation and rotation were 
performed with the same rejection criterioa For non-overlapping sites, the change in the potential 
energy accompanying the new configuration was calculated and subsequently accepted or rejected 
in accordance with standai-d acceptance probabilities (Molecular Simulation Inc., 1998). 
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Zeolite Modek 

The structure determined from powder neutron dijBfraction data by Fitch et aL 
(Fitch et al., 1986) was used as the structure for Na-Y Zeolite (Si/Al=2.43). The location of the 
Na cations were 32 atoms on Sn sites, 16 atoms on SI sites and 8 atoms on SI' sites. Na cations 
on SI and SF sites were located in such a manner that two adjacent SI and SI' were not 
simultaneously occupied. This is due to the electrostatic repulsion caused by their small separation 
of 0.218 nm (Fitch et al, 1986). The position of Al in Ihe framework was randomly chosen using 
the "disorder" function in the "crystal builder" module that obeyed the Loewenstein rule. No 
charge-free "blocking atom" was added because large molecules such as benzene were not 
located in the site, which was actually sterically inaccessible. 

For the simulation of Ag ion-exchanged zeolite (Si/Al=2.43), the Na cations in Na- 
Y were replaced by Ag maintaining the same cation locations. An all silica Y-Zeolite structure 
was used to simulate Ultra-Stable Y-Zeolite (Sj/A1=195) since this zeolite has less than one Al 
atom per unit cell. 

Until recent]^, the location of the Ag cation was not accurately determined. Using 
neutron diffraction data, Hutson et al (Hutson et al., 2000) reported the Ag location and 
occupancies of Ag in Ag-Y (Si/Al=2.43) after 450 °C dehydratioa Their results showed that the 
Ag cation sites were 28 atoms on Sn sites, 4 atoms on SII' site, 11 atoms on SI sites and 12 
atoms on SI' sites. Simultaneous occupancy of Ag cations at Sn and SII' as well as SI and SI' 
sites was unlikely due to the large repulsioa This location is very simiku* to the location of Na 
cations in Na-Y reported by Fitch et al (1986). Therefore, the zeolite model for Ag-Y used in this 
work is considered to reflect the actual Ag cation locations in Ag-Y (Si/Al=2.43). 
Forcefield Determination 

One of the most important fectors for GCMC simulation is the selection of an 
energy e?q>ression and parameters for the forcefield. The forcefield used in these simulations was a 
modified version of Cerius2 Sorption Dementis Forcefield, since Demontis et aL simulated the 
mobility of benzene in Na-Y zeolite successfiiUy using molecular dynamics methods (Demontis et 
al., 1989). In this model, the total potential energy (Uxotai) between the zeolite lattice and 
adsorbate molecules is written as the sum of the interactions between adsorbate molecules(AA) 
and that between the adsorbate molecules and the zeolite lattice (AZ). 
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Urotai^U^^U^^ (2) 
Interactions between adsorbate molecules (Uaa) were described as the sum of the atom-atom 
Bucldngham potentials (the first and second terms below), which have shown to yield reasonable 
predictions for liquid and solid benzene (Williams and Cox, 1984), and the Coulombic interaction 
potentials (the third term below). 



= 4 exp(-4.r^) — -I + 



Ml 



(3) 



The first and second terms represent, respectively, the repulsive and dispersiy^e potential energies 
between sites "i" and "j", and the third term represents the interaction potential between point 
charges qi and qj of sites "i" and "j" separated by a distance ry. 

The interactions between the adsorbate molecules and the zeolite lattice (Uaz) is 
written as the sum of the atom-atom Lennard- Jones 12-6 potential and the Coulombic interaction 
potentials. 
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Following the empirical approach used by Razman and Hall (Razmus and Hall, 
1991), Watanabe et al. (Watanabe et aL, 1995) used the empirical approach to parameterize the 
fprcefield for adsorption of N2, O2, and Ar in type A and fiiujasite zeolites. The dispersive part of 
the Leimard-Jones 12-6 potential (Udispeision) was written using an adjustable parameter (p.) of the 
adsorbate and the polarizabilities (a) of the atoms of the adsorbate molecules (i) and in the zeolite 
lattice 0)- 



The potential energy well depth parameter, Sy .is then written as 



r \ 



4<T 



(5) 



(6) 



The Lennard- Jones parameter, ay, for the adsorbate-zeolite interaction is gtvoi by the mixing rule, 
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The Oj for the zeolite lattice are related to the van der Waals radii, Bj, of the respective ions by 

/^=2^«o-^ (8) 
The well depth parameter, 8y , is assumed to be a geometric combination of 8] and bj, 
s,=(s,ejf' (9) 
The van der Waals energy in the periodic framework was calculated using a minimum image 
convention with an interaction cut-off distance of 1 nm. The Coulombic term was evaluated using 
the Evald summation method (Allen and Tildesley, 1987), 

The potential parameters used for Na-Y, Ag-Y and H-USY are listed in Table 15. 
The potential parameters for Na-Y are the same as that of Demontis et al (Demontis et al., 1989). 
Van d^ Waals interactions between adsorbate and zeolite cage were niodeled by means of 
interactions between the (C,H) and (Na,0) atoms. Interactions with the Si and Al atoms were 
neglected because th^ were well shielded by the o?grgen atoms of the Si04 and AIO4 tetrahedra. 
Electrostatic interactions between all charges were included in the potential energy calculation. As 
for Ag-Y and all silica Y zeolites, modification of parameters was performed using the mixing rule 
and the geometric combination rule. The detailed scheme for parameterization will be discussed 
shortly. 

Assignment of charges for adsorbate and adsorbent 

Assignment of charges also has a significant impact on simulation results when the 
Coulombic interactions are not negligible. In this work, the existing model for benzene in Cerius^ 
was used. In this model, each carbon in benzene carries -0.093 point charge, while hydrogen 
carries +0.093. These values were nearly the average values from the various charge 
determination methods: MOP AC (0.1022), ZENDO (0.072) and charge equilibration method (Qcq) 
by Rappe and Goddard (0.1147) (Rappe and Goddard, 1991). 

The point charges of zeolite were assigned either by (i) using the point charges in 
the literature or (ii) calculating by means of the charge equilibration method in the Cerius2 
software. In this work, Na, Si/Al, O in Na-Y carried the charges of +0.78, +1.2125 and -0.72, 
respectively, since these values were successfully used in quantum mechanics simulations 
(Kitagawa et al., 1996). Here, no distinction was made between Si and Al atoms in terms of point 
charges. Charges of Ag-Y and all-silica Y-Zeolites were calculated by the charge equilibration 
method: Ag(+0.45), Si/Ai(+1. 16875), O(-0.65), which will be referred to as Case I for Ag-Y, and 
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Si(+1.2650), O(-0.6325) for H-USY. In order to investigate the effect of the assigned charges to 
the isotherm simulation results, the same charges that were used for Na-Y, i,e,, Ag(-H).78), 
Si/Al(+1.2125) and O(-0.72), were abo used for Ag-Y. This case is referred to as Case H for 
Ag-Y. For all cases, regardless of the positions of cations and framework atoms (Si, Al, O), the 
same values of point charge were assigned to the same atoms throughout the simulation unit cell, 
even though the point charge might vary with the position of the atom. The Ag cation charges 
used for Monte Carlo simulation were also consistent with the values calculated by molecular 
orbital method at MP2/3-21G (Ag charge : +0.49) or MP2/LariL2DZ level (Ag charge : +0.74) 
using the zeolite clust^ model (Chen and Yang, 1996). 
Results and Discussion 
Experimental 

Benzene/cyclohexane adsorption isotherms 

Pure component isotherms of benzene and cyclohexane were measured for H- 
USY(Si/Al=195), Na.Y(Si/Al-2.43), Ag-Y(Si/Al=2.43), and Pd-Y(Si/Al=2.43) at 120 and 180 
^C, shown in Figures 22-25. Isosteric heats of adsorption calculated using the Clausius- 
Clapton equation are listed in Table 16. 

In the case of H-USY, the adsorbed amounts of cyclohexane were higher than that 
of benzene at > 0.04 atmosphere. This result qualitatively agreed with calculations using the 
Horvath-Kawazoe theory (Horvath and Kawazoe, 1983; Cheng and Yang, 1994). The threshold 
pore-filling vapor pressures of cyclohexane and benzene at 120 °C and 180 °C were calculated 
using the spherical pore model of the H-K method. The threshold pressure is the pressure where 
the steep rise in the isotherm occurs. At both temperatures, it was found that the pore-fiUing 
pressure of cyclohexane was 20 % lower than that of benzene. The reason for this difference was 
the larger polarizability (11.0x10"^"* cm^ for cyclohexane and 10. 3x10"^"* cm^ for benzene) and 
magnetic susceptibility (11.3x10"^^ cm^ for cyclohexane and 9.1x10"^^ cm^ for benzene) of 
cyclohexane. However, the calculated isosteric heat of adsorption of benzene was slightly higher 
than that of cycloh^cane. A possible reason for the larger heat of benzene adsorption was the 
^dstence of in the zeolite, although the number of proton was extremely small (less than one 
per xinit celQ. The interaction between protons in zeolite and 7c-electrons of benzene might be 
stronger than the interaction between oxygen in zeolite and benzrae. Therefore, the amounts of 
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benzene adsorption were slightly larger than that of cyclohexane at less than 0.04 atm, leading to 
its higher heat of adsorption. 

On all other sorbents investigated in this work, larger amounts of benzene were 
adsorbed at all pressure ranges. Heats of adsorption for benzene were 18 kcal/mol for Na-Y and 
21 kcal/mol for Ag-Y, while heats of adsorption for cyclohexane were 10-12 kcal/mol for botihi 
Na-Y and Ag-Y, These results are clearly due to the strong interactions betwem the cations in 
zeolite and the benzene molecule. Consequently, Ag-Y and Na-Y could adsorb large amounts of 
benzene even at the pressure of less than 0.01 atm. The heats of adsorption on Na-Y were 
comparable to previous work (Barthomeuf and Ife, 1973). Most importantly, Ag-Y showed the 
strongest interactions for benzene and the largest heats of adsorption (20.5 kcal/mol). For 
cyclohexane adsorption, the adsorption amounts and heat of adsorption were almost the same 
between Ag-Y and Na-Y, since Tc-complexation was absent for cyclohexane. A detailed discussion 
will be made in the section of molecular orbital calculation and Monte Carlo simulation results. 
These results indicated that improved purification performance of benzene from cyclohexane is 
possible using Ag-Y. The isotherms of both Ag-Y and Na-Y are reversible, although it took a 
little longer time to desorb benzene and cyclohexane. For example, 200 minutes are necessary to 
desorb benzene completely after the adsorption at 0.1 atmosphere. This slow desorption is clearly 
due to the strong interaction between benzene and Ag, However, the bonds are of the order of 20 
kcal/mol, which are weak enough so a substantial fraction of the adsorbed benzene is desorbed in 
a few minutes. For practical applications, thermal desorption may be a desirable option. 

Pd-Y gave lower adsorption amounts than both Ag-Y and Na-Y. This decrease of 
adsorption amounts by ion-exchange with Pd was due to degradation of the zeolite structure. In 
Figure 26, the pore size distributions for all sorbents calculated from N2 isotherms at 77 K are 
shown. Pd-Y had a small pore volume with a wide pore size distribution, clearly indicating the 
collapse of the zeolite pore structure or pore blockage by Pd. The Pd ion-exchanged was 
performed at 80 ^C, compared with the Ag ion-exchange at room temperature. Collapsing of 
zeolite structure tends to occur considerably more in acidic solution and at higher temperatures. 
The pore volimie of Ag-Y also decreased by 23% after ion-exchange. However, this reduction of 
Ag-Y was mainly caused by 37% density increase of the zeolite lattice by replacing Na"^ with Ag^ 
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(27% reduction of pore volume per weight basis). The pore structure of Ag-Y was intact after 
ion-exchange. 

Another interesting phenomenon in Figure 26 is the apparent shift of the pore size 
distribution in H-US Y toward a higher value, compared to other zeolites. This was likely due to 
the lowK- polarizability of oxygen in the framework of high-Si02 zeolites (Pellenq and Nicholson, 
1993), which will be discussed in detail in the section of Monte Carlo simulation results. 
Otherwise, the interactions between nitrogen and cations in zeolites would have led to a shift of 
the pore sizes toward smaller values for Na-Y, Ag-Y and Pd-Y. Since the pore sizes of Na-Y, 
Ag-Y and Pd-Y were the same as the actual value for Y zeolites (1.3nm), the latter explanation is 
not plausible. 
Separation factors 

The pure component adsorption data were first fitted by both D-A (Dubinin- 
Astakhov) and L-F (Langmuir-Freundlich) isotherms. The separation fectors of benzene over 
qrclohexane were then obtained using the equivalent multi-component isotherms. In the D-A 
equation of pure component gas, the volume adsorbed, V, at relative pressures P/Ps is expressed 
as 

^^=^iexp[-[cin^) (10) 

where 

C = — (11) 

Doong and Yang (Doong and Yang, 1988) extended the D-A equation to mixed-gas adsorption in 
a simple way by using the concept of maximum available pore volume without any additional 
equations such as the Lewis relationsh^ (Lewis et aL, 1950). For binary mixtures. 




f^=(F„-F,)exp 



^'2=(^o2-I^)expl-lQln^ 



(12) 



(13) 



Since equation (12) and (13) are linear, Vi and V2 can be obtained in a straightforward fiishion. 
After the calculation of Vi and V2, volumetric adsorbed amounts were converted to molar 
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adsorbed amounts and the separation fectors were calculated. The conversion of molar adsorbed 
amount to volumetric adsorbed amount and vice versa was performed using the liquid density 
data of the adsorbates. The liquid density at 120 and 180 were obtained by the modified 
Radcett equation by Spencer and Adler (Spencer and Adler, 1978). The liquid densities used in 
this work w^e 0.769 g/cc at 120 °C and 0.689 g/cc at 180 ° C for benzene, and 0.679 g/cc at 120 
°C and 0.606 g/cc at 180 for cyclohexane. The adsorbed molar volumes were assumed to be 
the same in both single- and mixed-gas adsorptioa 

For the Langmuir-Freimdlich isotherm equation, pure component isotherm is 
written using three adjustable parameters: 

TheL^i'' equation can also be extended for n-component mixture (Yang, 1987): 

<,.-^4^ (15) 

The fitting parameters for both isotherms are listed in Table 17 and the fitted 
curves are plotted in Figures 22 - 25. The separation fectors of benzene over cyclohexane were 
calculated using the equilibrium mole firactions in adsorbed and gas phases: 

rv = "^fegw^gwg / ^benzene 

X lY ^ ^ 

cyclohaaaie f cydohexane 

The calculated separation factors for Ag-Y and Na-Y are compared in Table 18. Both D-A and L- 
F equations led to similar dependence of separation factor on pressure. The separation &ctor for 
Ag-Y increased with decreasing concentration of benzene, while that for Na-Y at 180 °C 
decreased. Although the separation Motors calculated by the two isotherms were not the same, 
the advantages of Ag-Y over Na-Y are clearly seen, especially at lower concentrations of 
benzene. This is certainly because Ag-Y can maintain its adsorption ability to benzene even at 
aroimd 1x10'^ atm. On the other hand, at 0.01 and 0.001 atm of benzene and at the lower 
temperature of 120**C, Na-Y showed equivalent separation fectors to Ag-Y. 
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For separations by adsorption and membranes, a separation fiictor of the order of 
10 is considered high. For example, in air separation by pressure swing adsorption using zeolites 
and also by polymeric membranes, systems with separation fiictors (N2/O2 for adsorption and 
O2/N2 for membranes) of near 10 are commercially used. For aromatics/aliphatics separations 
iising liquid membranes (e.g., toluene/heptane and benzene/heptane separations), separation 
fectors of the ord^ of 10-100 yielded excellent separations (Li, 1968; Li, 1971). The separation 
&ctors shown in Table 18, in the order of 10-10,000, are excellent for separation/purification. 

Diffusion time constants 

In Table 19, diflftision time constants calculated from 0-0.3 fractional uptake were 
summarized. Diffusion time constants of benzene and cyclohexane in Na-Y and Ag-Y were of the 
order of 10"^ 1/s. These values were almost the same as those of benzene and cyclohexane into 
AgNOs monolayer dispersed Si02 sorbent (Takahashi et aL). However, these values were an 
order of magnitude lower tiian those of propylene and propane in AgNOa monolayer-dispersed 
Si02 sorbent (Rege, 1998) and at least the same as those of ethylene and ethane into Ag ion- 
exchanged resins (Wu et aL, 1997). The difiijsion rates of benzene in these zeolites are high for 
applications (Yang, 1987). 
Molecular Orbital Calculation 
Adsorption Bond Energy and Molecular Geometries 

Adsorption energy and structural geometry were calculated using molecular orbital 
theory. The calculated heat of adsorption using B3LYP/3-21G level by equation (1) is also 
compared with experimental and Monte Carlo simulation results in Table 16. The theoretical 
value, 17.2 kcal/mol, was in fiiirly good agreement wdth the experimental value. The comparison 
of bond length by benzene adsorption revealed that the Ag-O length inoreased from 3. 023 A to 
3.246A. The C-C bond length in benzene also increased from 1,385A to 1.385-1,391 A. This trend 
was consistent with that for benzene interaction with AgCl (Takahashi et al,). 
Natural Bond Orbital (NBO) Results 

The nature of benzene-Ag zeolite interaction can be better understood by 
analyzing the NBO results. The main feature of the interaction can be seen from the population 
changes in the vacant outer-shell s orbital of the silver and those in the d orbitals of silver upon 
adsorption. In Table 20, NBO analysis results for benzene on Ag zeolite are compared with other 
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combinations of adsorbates and adsorbents. According to the electron population changes by 
benzene adsorption on Ag zeolite, the trend is in agreement with previously reported trend of tc- 
comple?cation (Cheng, 1996; Huang et al, 1999a; 1999b; 1999c ; Takahashi et al.) i.e., donation 
of electron charges from the p-orbitals of benzene to the vacant s-orbital of the silver (a donation) 
and, simultaneously, back donation of electron charges from the d-orbitals of silver to tc* orbital 
of benzene (d-7C* back donation). The contribution of d-7c* back donation is larger than a 
donation, resulting in the decrease of net electron charge in silver. The net electron charge of 
carbon in benzene increased instead. 

The differences of electron population changes between benzene and ethylene 
adsorption deserve a discussioa From Table 20, it is seen that the contribution of a donation is 
predominant in ethylene adsorption, while d-TU* back donation is slightly larger in benzene 
adsorptioa Also, the extent of a donation and d-Tc"^ back donation by benzene adsorption was 
considerably smaller compared with ethylene adsorptioa Therefore, 7c-complexation with benzene 
was weaker than that with ethylene. Weak Tc-complexation of benzene can also be understood by 
the difference in bond lengths. The bond length between carbon in ethylene and silver in AgCl was 
2.66 A (Huang et al., 1999b), while that between carbon in benzene and silver in Ag zeolite or 
AgCl was 2,71-3.7 A or 3.22 A (Takahashi et al.), respectively. This longer bond length was 
caused by the less overlap of orbitals from the weaker 7i-complexation. 
Monte Cario Simulation 

Monte Carlo simulation of adsorption isotherms provides great insight into the 
adsorption mechanism via model construction, charge assignment and forcefield determination In 
this section, benzene adsorption on Na-Y, H-USY and Ag-Y is simulated and the mechanism of 
adsorption on each zeolite is discussed. 
Na-Yzeolite 

GCMC simulation results of benzene isotherms at 120 and 180^C using the 
potential parameters given in Table 15 are shown in Figure 27, where they are also compared with 
experimental data. The simulation results were in excellent agreement with experimental data. The 
calculated heats of adsorption were 16.4 - 18.5 kcal/mol, which were in agreement with that from 
experiments. It was hence confirmed that the energy expressions and parameters were suitable for 
simulation of benzene isotherms on Na-Y. 
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H'USYzeolite 

First, simulation of benzene isotherms on H-US Y was att^pted using the same 
potential parameters as on Na-Y. However, the simulated isotherms were considerably higher 
than the experimental data. This disagreement was caused by incorrect well depth parameters, 8cm: 
and so-H, for van der Waals interaction. Pellenq and Nicholson (Pellenq and Nicholson, 1993) 
reported the values of the polarizabilities of firamework aton:is in silicalite and ahiminosilicates 
detennined from Auger electron spectroscopy data. It was found that the polarizabihty for oxygen 
in zeoite framework is sensitive to the ratio of silicon to aluminum and also the nature of the 
extraframework cations. Watanabe et al (Watanabe et al, 1995) used different potential values 
for oxygen in silicalite (0.202 kcal/mol), denoted by Oz(Si-O-Si), and that in aluminosilicates 
(0.334 kcal/mol), denoted by Oz(Si-0-Al). 

In this work, an adjustmrat of soz(si-o^i)-c and eQz(si-o^i>H was made to simulate the 
experimental Henry's law constants at 120 and 180 ®C. Here since low pressure data depend only 
on the benzene-zeolite interactions so the benzene-benzene interactions were neglected. Using 
the well depth parameter for 8Na^ = 0.04187 kcal/mol and 8Na-H= 0.03153 kcal/mol for Na-Y (in 
Table 15) and assuming ENa= 0.041 kcal/mol as used hy Watanabe et al. (Watanabe et al, 1995) 
for Na-A and Na-X, sc = 0.0428 kcal/mol and Eh = 0.0243 kcal/mol were obtained from the 
geometric combination rule, equation (9). Therefore, SQz(si.o-si)-c and soz(si.o-si)-H for H-USY could 
be expressed as (0.0428 x 8cte(si-o-si))^'^ and as (0.0243 x soz(fflo-si))^*^ , respectively. By fitting with 
ecte(si<«i), eQz(si-o-si) = 0.665 kcal/mol was obtained to fit the experimental Hrary's law constant of 
0.519 molecules/(cell.kPa) at 120 and 0.137 molecules/(cellkPa) at 180 The ratio of SQz<a- 
osi) over 8Qz(si^Ai) was 0,44 (calculated from BQz(sio-Ai)-H ) - 0.60 (calculated from 8Qz(si-o-Ai)-c ), 
whidi was in agreement with the ratio of 0.60 obtained by Watanabe et at (Watanabe et al., 
1995). 

Using the modified parameters for H-USY (in Table 15), the benzene isotherms 
were calculated and compared with experimental isotherms in Figure 28. Again, the isotherms and 
the heat of adsorption (8.2 - 9.2 kcal/mol) obtained from simulation were in excellent agreement 
with the experimental data. 
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Ag-Y zeolite 

The potential parameters for Ag-Y are listed in Table 15. These parameters were 
modified to get the low pressure points of the isotherm at 180 °C (1.30 mmol/g at 1.6 x 10"^ 
atmosphere), although adjustment of the parameters to fit the experimental Henry's law constant 
was preferable. The reason not to use Ihe e?qperimental Henry's constant was that the small 
pressure values ( 8.7 x 10"^ atmosphere) at the lowest adsorbed amount in the isotherms might 
lead to the large error of potential parameters. The procedure of adjustment was basically the , 
same as in the case for H-USY except for changing the SAg value instead of the Soz(si-o-si) value. Sc 
= 0.0428 kcal/mol and 8h = 0.0243 kcal/mol were obtained from the geometric combination rule, 
equation (9). Then, SAg^ and eAg-nfor Ag-Y could be estimated from (0.0428 x 8Ag)*^^ and (0.0243 
X eAg)°*^, respectively. In Case I for Ag-Y, SAg became 14.3 kcal/mol, while eAg= 13.5 kcal/mol 
was obtained for Case DL The larger well-depth parameter was needed for Case I to compensate 
for the weaker Coulombic interaction due to the smaller Ag cation charge. In both cases, the &Ag 
values used were more than 2 orders of magnitude greater than SNa, clearly indicating the strong 
interaction from Tc-complexatioa 

Figure 29 shows the simulation results of benzene on Ag-Y for Case I and Case n. 
The differences between Case I and Case n were feirly small. In both cases, the agreement with 
experimental isotherms was not as good as that for Na-Y and H-USY. However, they certainly 
followed the basic trend of the experimental isotherms. The calculated heat of adsorption was 
18.8 - 21.1 kcal/mol, which was higher than Na-Y and also in agreement with the experimental 
value. 

In all cases, not only Ag-Y but also Na-Y and H-USY, the simulated isotherms 
traded to be higher than the experimental isotherms in the lower pressure range with few 
exceptions, while they were lower in the higher pressure range. The higher simulation results at 
low pressures might be caused by the larger values for the well-depth parameter between oxygen 
(zeolite) - carbon (benzene) interaction and oxygen (zeolite) - hydrogen (benzene) interaction. 
The other possible cause was the imperfection of the zeolite crystals. Since zeolites contain 
defects vs4iile in Monte Carlo simulation perfect zeolite lattice was assumed. In the high pressure 
range, larger experimental values were possibly caused by adsorption between zeolite crystals. 
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Adsorption sites for benzene on Na-Y and Ag- Y were investigated using the Mass 
Cloud Plots function in Cerius2 (not shown here). It was confirmed that benzene was adsorbed 
preferentially near the SII cation sites and the center of the 12-ring window, which was consistent 
with the neutron difl&action data for Na-Y obtained by Fitch (Fitch et aL, 1986). For Ag-Y, 
benzaie was adsorbed slightly further away from the Sn cation sites than that in Na-Y, due to thie 
feet that the van der Waals radius of Ag* is 30% larger than that of Na*. 
Notation 

A = repulsive parameter for Buckingham potential 

B = repulsive parameter for Buckingham potential or Langmuir constant 

C = dispersive parameter for Buckingham potential or contant for D- A equation 

D = difiiisivity 

E = energy 

n = constant 

P = pressure 

q - point charge or molar adsorbed amount 
R = van der Waals radii or gas constant 
r = distance or radius 
T = temperature 
U = potential energy 

V = volumetric adsorbed amount 

X = equilibrium mole fraction in adsorbed phase 

Y = equilibrium mole fraction in gas phase 
Greek letters 

a = polarizability or separation factor 
P = adjustable parameter 

8 = well depth paramet^ for Lennard- Jones potential 
a = distance parameter for Lennard- Jones potential 
Subscript 

i = atom site or componet 
j = atom site or component 
s = saturation 
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Table l5The Potential Parameters used for Na-Y, Ag-Y and H-USY 



Na-Y 






CO 


H-H 


C-H 








A (Kcal/mol) 


88371 


28ol 










D /I /AN 




^ 7 A 

o,/*f 


o.o / 








C (A^*kcal/mol) 


583.13 


32.60 


137,88 












vj-n 


Wa.P 


Ma.M 






8 (Kcal/mol) 


0.25547 


0.16515 


0.04187 


f\ rtO ICO 

0.03153 








2.9984 


2.5893 


3.3279 


2.7651 


AgV 






0-C 


0-H 


Ag-C 


Ag-H 


(Case 1) 


1 


B (kcal/mol) 


0.25547 


0.16515 


0.7823 


0.5895 






crA) 


2.9984 


2.5893 


3.5773 


3.0745 


Ag-Y 


Uaz 




0-C 


0-H 


Ag-C 


Ag-H 


(Case II) 




e (kcal/mol) 


0.25547 


0.16515 


0.7607 


0.5728 








2.9984 


2.5893 


3,5773 


3.0145 


H-USY 


Uaz 




0-0 


O-H 







8 (kcal/mol) 0.1687 0.1271 

siA} 2.9984 2.5893 

1) Parameters for U aa are the same for all sorbents. 

2) The numbers with underline were determined in this work. 
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Table Comparison of Heat of Adsorption (kcal/mol) 
between Experimental and Simulated Results 



Na-Y Ag-Y H-USY 



Benzene Experimental 17.7-18.2 20.5-20.7 6.2-7.2 

Monte Cairo Simulation 16.4-18.5 18.8-21.1 8.2-9.2 
Molecular Orbital Calculation *** 17.2 

• Barthomeuf et al. ^> 17.518.0 *** 

Cyclohexane Experimental 10.9-13.0 12.1-12.3 4.1-4.3 

Barthomeuf et al. ^> 12.0-13.2 *** 



l)Barthomeuf and Ha, 1973 
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Table|7 Parameters for the D-A and L-F Equations 











D-A Equation 






L'r iiC|uaiion 




AdsorDent 


Adsorbate 


1 emp. V W 


*o 


0 


n 


PQ 




B 


n 












CO- 


^atm) 


inrimoi/Mj 


I X/a1.lTl J 


t "1 


IMa-Y 


Benzene 


1 on 


n ooB% 


0.092 


3.76 




1 


30.7 


2.53 




Benzene 


1 on 
loU 




0.120 


4.73 


1 n 07 




t'l ^ 


0.9 




Cyclohexane 


1 on 


n ofi7 


0.149 


5.01 




^ . o o 


51.2 


1.18 




f^f/^l rsh Avon A 

v^yciQneAane 


X ov 




0.168 


6.44 


9.48 


2.30 


61.9 


0.77 


Aa V 


Denzene 


1 9n 




0.056 


1,57 


2.96 


3.07 


3.67 


7.46 




Benzene 


loU 




0.082 


3.05 


1 n cn 

XU.U / 






1 85 




Cyclohexane 


120 


0.225 


0.133 


4.32 


2.84 


2.02 


48.4 


1.44 




Cyclohexane 


180 


0.237 


0.145 


4.96 


9.48 


1.80 


74.3 


0.91 


H-USY 


Benzene 


120 


0.237 


0.318 


2.84 


2.96 


3.40 


11.2 


0.62 




Benzene 


180 


0.234 


0.333 


2.17 


10.07 


3.40 


1.11 


0.74 




Cyclohexane 


120 


0.249 


0.265 


5.85 


2.84 


1.33 


2500000 


0.2 




Cyclohexane 


180 


0.15 


0.232 


4.93 


9.48 


0,73 


260 


0.39 


Pd-Y 


Benzene 


120 


0.173 


0.237 


2.15 


2.96 


0.883 


240 


0.7 




Cyclohexane 


120 


0.16 


0.317 


1.34 


2.84 


0.753 


9.22 


1.2 



58 



\VO 03/020850 PCT/US02/28040 



Table Separation Factors estimated by the D-A and L&F Equations 

IZ 



Pressure (atm) 




D-A Equation 


L-F Equation 


Benzene Cvclohexane 


TenriD. 


ar-y _ 


In3-Y 


Arr V 

Ag-Y 


Na Y 


0.01 1 


120 


107 




6 


12 


0.001 1 


120 


416 


441 


45 


48 


0.0001 1 


120 


2286 


1016 


335 


191 


0.00001 1 


120 


14361 


2553 


2462 


767 


0.01 1 


180 


39 


95 


19 


5.4 


0.001 1 


180 


120 


71 


57 


4.2 


0.0001 1 


180 


456 


26 


163 


3.3 


0.00001 1 


180 


1767 


0.5 


470 


2.6 
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Table Diffusion Time Constants of Benzene and Cydohexane 



Adsorbents 


emp. 


Benzene 


(Pressure Change) 


Cyclohocane 


(Pressure Change) 




( C) 


(1/s) 


(atm) 


fl/s) 


(atm) 


Na-Y 


120 


2.1x10-* 


(5.4x1 0-'->2.6x10-^ 


Q ftvin*^ 

9.0X 1 \f 


a 1x10"^— >2 0x10"^ 


Ag-Y 


120 


4.2x10-* 


(5.4x1 0-*->2.6x10-^ 


1.1x10-* 


(1.1x10-'->2.0x10-*) 


AgNQ3/SiO2(0.33g/g) 


120 


1.3x10-* 


(2.0x1 0-'->3.7x10-^ 


2.9x10-* 


(2.0x1 0-«->3.7x10-^ 


Adsorbents 


Temp. 


ropylene (Pressure Change) 


Propane 


(Pressure Change) 




(-C) 


(1/s) 


(atm) 


(Ife) 


(atm) 


AgNQ3/SiQ2(0.33g/g) 


70 


3.5x10-* 


(0->1) 


*** 




Ag-Amberlyst-35 ®' 


70 


LSxIO-* 




1.4x10-^ 


*** 



1) Takahashietal..2000 

2) Regeetal..1998 

3) Wuetal., 1997 
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Table Summary of the NBO Analysis of then-complexation: 
Electron Population Changes of Each Orbital after Adsorption 









Silver 






Carbon 






Method / BasisSet 


5S 


S4d 


Net Chance 


2s 


£2p 


etChanae 


enzene on Ag Zeofite 


B3LYP/3-21G 


0.00541 


-0.0135 


-0.00805 


0.00554 


0.00709 


0.01263 


Benzene on AgCI ^ 


B3LYP/LanL2DZ 


0.0033 


-0.0071 


-0.0038 








C2H4 on Ag Zeolite ^ 


WIP2/3-21G 


0.0596 


-0.0177 


0.0419 


0.0203 


0.0193 


0.0396 


C2H4 on AgCI ^ 


MP2/3-21G 


0.1204 


-0.0470 


0.0734 


0.0226 


0.0142 


0.0368 


C2H4on AgCI^ 


B3LYP/LanL2DZ 


0.061 


-0,055 


0.006 


*** 


*** 





1) Takahashietal..2000 

2) Chen and Yang, 1996 

3) Huang etal., 1999 
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KXAMPLES OF SOME SUI T ART F. ADSORBENTS 

Abrief description of some non-limitative examples of adsorbents 
which may successfully be used in the present invention follows. Detailed 
descriptions may be found in U.S. Patent No. 6,423,881, and in U.S. Patent No. 
6,215,037, each of which applications is incorporated herein by reference in its 
entirety. 

The adsorbent comprises a carrier having a surface area, the carrier 
having dispersed thereon a monolayer of a d-block transition metal compound, eg. a 
silver compound. The silver compound releasably retains the thiophenes; and the 
carrier comprises a plurality of pores having a pore size greater than the eflfective 
molecular diameter of the thiophenes. 

It is to be understood that any suitable carrier may be used. However, 
in a prefeired embodiment, the carrier has a BET surfece area greater than about 50 
square meters per gram and up to about 2,000 square meters per gram, and comprises 
aiplurality of poi«s having a pore size greater than about 3 angstroms and up to about 
10 microns. In a more preferred embodiment, the carrier is a high surfece area 
support selected ftom the group consisting of refractory inorganic oxide, molecular 
sieve, activated carbon, and mixtures thereof Still more preferred, the carrier is a 
refractory inorganic oxide selected from the group consisting of pillared clay, 

damina and silica. 

It is also to be understood that any suitable silver compound may be 
used. However, in a preferred embodiment, the silver compound is a silver (T) 
halide. In a more preferred embodiment, the silver compound is a silver salt, and the 
salt is selected from the group consisting of acetate, benzoate, bromate, chlorate, 
lasrchlorate, chlorite, citrate, fluoride, nitrate, nitrite, sulfete, and mixtures thereof 

In one exemplary embodiment of this embodiment of the present 
invention, the silver compound is silver nitrate (AgNOs) and the carrier is silica 
(SiOz). 

The method of the present invention may further comprise the step of 
dtan^ng at least one of the pressure and temperature to thereby release the 
thiophenes-rich component &om the adsorbent. It is to be understood that the 
pressures and temperatures used may be witiiin a suitable range. However, in the 
preferred embodiment, the selected pressure of preferential adsorption is a first 
pressure, and the pressure of release is a second pressure less than the first pressure. 
£ea more preferred embodiment, the first pressure is in a range of about 1 
atmosphere to about 35 atmospheres, and the second pressure is in a range of about 
0.01 atm to about 5 atm. 
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In the preferred embodiment, the selected temperature of preferential 
adsorption is a first temperature, and the temperature of release is a second 
temperature greater than the first temperature. In a more preferred enibodiment, the 
first temperature is in a range of about CC to about SCC, and the second 
tesnperature is in a range of about 70°C to about 200 ^'C. 

Without being bound to any theory, it is believed that the retaining of 
the thiophenes is accomplished by formation of k -complexation bonds between the 
silver compound and the thiophenes. Separation by n -complexation is a subgroup of 
chemical complexation where the mixture is contacted with a second phase, which 
contains a complexing agent. The advantage of chemical complexation is that the 
bonds formed are stronger than those by van der Waals forces alone, so it is possible 
to achieve high selectivity and high capacity for the component to be bound. At the 
same time, the bonds are still weak enough to be broken by using simple engineering 
operations such as raising the temperature or decreasing the pressure. 
15 The 7c -complexation generally pertains to the main group (or d-block) 

triansition metals, that is, fi-om Sc to Cu, Y to Ag, and La to Au in the periodic table. 
These metals or their ions can form the normal o bond to carbon and, in addition, the 
unique characteristics of the d orbitals in these metals or ions can form bonds with 
the unsaturated hydrocarbons (olefins) in a nonclassic manner. This type of bonding 
isd)roadly referred to as tc -complexation, and has been considered for gaseous 
hydrocarbon separation and purification using cumbersome liquid solutions. 

A fiirther example of an adsorbent comprises an ion-exchanged 
zeolite selected fi:'om the group consisting of zeolite X, zeoUte Y, zeolite LSX, and 
mixtures thereof, the zeolite having exchangeable cationic sites, and at least some of 
ti® sites have a d-block transition metal cation, eg. silver cation or copper cation, 
present. In the more preferred embodiment, the cationic sites of the ion-exchanged 
zeolite contain silver cation. In a further preferred embodiment, essentially all 
cationic sites of the ion-exchanged zeolite contain the silver cation. 

Without being bound to any theory, it is believed that the preferential 
^tteorption occurs by 7i-complexation. 

Figs. 30-32 are graphs depicting fiirther breakthrough curves. 

While preferred embodiments of the invention have been described in 
detail, it will be apparent to those skilled in the art that the disclosed embodiments 
may be modified. Therefore, the foregoing description is to be considered exemplary 
lusher than limiting, and the true scope of the invention is that defined in the 
following claims. 



63 



wo 03/020850 



PCT/US02/28040 



What is claimed is: 

1 . A process for removing thiophene and thiophene compounds from liquid 
fiiel, the method comprising the step of: 

contacting the liquid fixel with an adsorbent which preferentially adsorbs the 
thiophene and thiophene compounds, at a selected temperature and pressure, thereby 
producing a non-adsorbed component and a thiophene/thiophene compound-rich 
adsoibed component, wherein the adsorbent comprises an ion-exchanged zeolite 
selected from the group consisting of zeolite X, zeolite Y, zeolite LSX, and mixtures 
thereoi^ tiie zeolite having exchangeable cationic sites, and at least some of the sites 
having d-block transition metal cation present, and wherein the preferential 
adsorption occurs by Tc-complexation. 

2. A process for removing thiophene and thiophene compounds from liquid 
£usbl, the method comprising the steps of: 

contacting the liquid fiiel with an adsorbent which preferentially adsorbs the 
thiophene and tiiiophene compounds, at a selected temperature and pressure, thereby 
producing a non-adsorbed component and a thiophene/thiophene compound-rich 
adsorbed component; the adsorbent comprising a carrier having a surface area, the 
(20Tier having a monolayer of a d-block transition metal compound dispersed on 
substantially the entire surface area, the metal compound releasably retaining the 
thiophene/thiophene compounds; and the carrier comprising a plurality of pores 
having a pore size greater than the effective molecular diameter of the 
thiophene/thiophene compounds; and 

25 changing at least one of the pressure and temperature to thereby release the 
thiophene/thiophene compound-rich component from the adsori>ent. 

3. The process as defined in claim 1 wherein at 10"^ atm, the adsorbent 
adsori^s more than 1 mmol/gram of thiophene. 

30 

4. The process as defined in claim 2 wherein at 10'^ atm, the adsorbent 
adsorbs more than 1 mmol/gram of thiophene. 

5. A process for removing aromatic compounds from liquid fiiel, the method 
caanprising the step of 

contacting the liquid fiael with an adsorbent which preferentially adsorbs the 
aromatic compounds, at a selected temperature and pressure, thereby producing a 
non-adsorbed component and an aromatic compound-rich adsorbed component. 
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wherein the adsorbent comprises an ion-exchanged zeolite selected from the group 
consisting of zeolite X, zeolite Y, zeolite LSX, and mixtures thereof, the zeolite 
having exchangeable cationic sites, and at least some of the sites having d-block 
transition metal cation present, and wherein the preferential adsorption occurs by ti- 
c^plexation. 

6. A process for removing aromatic compounds from liquid fuel, the method 
comprising the steps of: 

10 contacting the liquid fiiel with an adsorbent which preferentially adsorbs the 
aromatic compounds, at a selected temperature and pressure, thereby producing a 
non-adsorbed component and an aromatic compound-rich adsorbed component; the 
adsorbent comprising a carrier having a surface area, the carrier having a nipnolayer 
of a d-block transition metal compound dispersed on substantially the entire surface 
ats^a, the metal compound releasably retaining the aromatic compounds; and the 
carrier comprising a plurality of pores having a pore size greater than the effective 
molecular diameter of the aromatic compounds; and 

changing at least one of the pressure and temperature to thereby release the 
aromatic compound-rich component from the adsorbent. 

20 
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